AD- A 189  237 


m  FILE  COP" 


SACLANTCEN  REPORT 
serial  no.:  SR-112 


SACLANT  ASW 
RESEARCH  CENTRE 
REPORT 


t [  •  if' 


DTIC 

*ELECTE 

.  DEC  1  1 1967 


Low-frequency  reverberation 
measurements  with  an  activated 
towed  array: 

Scattering 
strengths  and  statistics 


D.  Marandino 


October  1987 


The  SACLANT  ASW  Research  Centre  provide*  the  Supreme  Allied 
Commander  Atlantic  (SACLANT)  with  scientific  and  technical 
aa »i»t an ci  under  the  term*  or  It*  NATO  charter,  which  entered  into 
force  on  1  February  1903.  Without  praludke  to  this  main  task — 
and  undar  the  policy  direction  of  SACLANT— -the  Centre  also  renders 
scientific  and  technical  assistance  to  the  Individual  NATO  nations. 


Approved  far  pobtte  rein 
Dtstrlbudoc  UnSmHad 


This  document  is  released  to  a  NATO  Government 
at  the  direction  of  SACLANT  ASW  Research  Centre 
subject  to  the  following  conditions: 

•  The  recipient  NATO  Government  agrees  to 
use  its  best  endeavours  to  ensure  that  the 
information  herein  disclosed,  whether  or  not  it 
bears  a  security  classification,  is  not  dealt  with 
in  any  manner  (a)  contrary  to  the  intent  of 
the  provisions  of  the  Charter  of  the  Centre, 
or  (b)  prejudicial  to  the  rights  of  the  owner 
thereof  to  obtain  patent,  copyright,  or  other  like 
statutory  protection  therefor. 

•  If  the  technical  information  was  originally 
released  to  the  Centre  by  a  NATO  Government 
subject  to  restrictions  clearly  marked  on  this 
document  the  recipient  NATO  Government 
agrees  to  use  its  best  endeavours  to  abide  by 
the  terms  of  the  restrictions  so  imposed  by  the 
releasing  Government. 


Page  count  for  SR-112 
(excluding  covers) 


Pages  Total 

i-iv  4 

1-44  44 

@1~@38  38 

86 


SAC’LANT  ASW  Research  Centre 
Viale  San  Bartolomeo  400 
19026  San  Bartolomeo  (SP),  Italy 

tel:  0187  540  111 
telex.  271148  SACENT  1 


WORTH  ATLANTIC  TREATY  ORGANIZATION 


SaCLANTCEN  SR-112 


Low-frequency  reverberation 
measurements  with  an 
activated  towed  array: 

Scattering  strengths 
and  statistics 


D.  Marandino 


The  content  of  this  document  pertains 
to  work  performed  under  Project  02  of 
the  SACLANTCEN  Programme  of  Work. 
The  document  has  been  approved  for 
release  by  The  Director,  SACLANTCEN 


Ralph  R.  Goodman 
Director 


*  Original  contains  color  \ 
plates:  All  DTIC  reproduet- 
I 'ns  will  be  In  black  and 

*>  te* 


Aooesslon  For 


HTIS  GRAJtl 
DTIC  TAB 
Unannounced 
Justification- 


D 

□ 


Distribution/ 


Availability  Codes 
[Avail  and/or 
Special 


Diet 


Saclantcen  SR 


intentionally  blank  page 


SACLANTCEN  SR- 11 2 


Low-frequency  reverberation  measurements 
with  an  activated  towed  array: 

Scattering  strengths  and  statistics 
D.  Marandino 


A 

Abstract:  Results  of  reverberation  measurements  made  with  a  low 
frequency  activated  towed  array  sonar  system  in  two  deep-water  locations 
of  the  Mediterranean  Sea  are  presented.  Measurements  took  place  under 
a  variety  of  environmental  conditions:  the  frequencies  used,  in  the  sub 
kiloherts  regions,  were  in  two  bands  around  370  and  740  Hr  Monostatic 
reverberation  data  were  collected  with  a  variety  of  waveforms  (pulsed  con 
tinuous  wave  and  linear  frequency  modulated  pulses)  in  order  to  obtain  an 
adequate  sample  of  the  scattering  function.  The  build-up  of  the  data  base 
has  stressed  data  collection  under  possibly  typical  operational  condition* 
for  further  performance  analysis  studies. 

Preliminary  findings  of  the  data  analysis  indicate  that  volume  reverberation 
does  not  seem  to  be  significant  for  this  experimental  set-up  and  that  surface 
reverberation  appears  almost  always  overshadowed  by  seafloor  scattering. 
Bottom  scattering  dominates  returns  through  both  the  so-called  fathometer 
effect  and  direct  backscattering  from  the  range  of  grating  angles.  Scatte¬ 
ring  strengths  appear  independent  of  frequency  and  spectral  histories  show 
no  significant  doppler  shift.  Long-range  returns  are  dominated  by  bottom 
highlights.  Reverberation  time  series,  which  are  highly  non-stationary,  can 
however  be  assumed  to  be  locally  covariance  stationary,  and  logarithmic 
normalisation  techniques  could  usefully  be  applied.  ^ 
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1.  Introduction 

This  report  presents  results  of  reverberation  measurements  made  with  an  activated 
towed  array  and  a  related  analysis  based  on  a  comparison  with  model  predictions. 
The  measurements  were  made  in  various  deep  water  locations  of  the  Mediterranean 
Sea  and  under  a  variety  of  sea-states.  The  transmisssion  frequencies  used,  in  the 
subkilohertz  region,  were  located  around  two  bandwidths  centered  at  370  Hz  and 
740  Hz.  Monostatic  reverberation  data  were  collected  using  a  variety  of  waveforms 
(continuous  wave  and  frequency  modulation  types)  in  order  to  obtain  an  adequate 
sample  of  the  reverberation  scattering  function. 

The  task  was  undertaken  in  support  of  on-going  studies  at  SACLANTCEN  which 
are  investigating,  both  theoretically  and  experimentally,  the  potential  for  long-range 
detection  using  a  sonar  based  on  the  low-frequency  activated  towed  array  concept. 
As  with  any  high-power  active  sonar  system,  reverberation  represents  a  serious 
system  performance  limitation  which  must  be  taken  into  account. 

Consequently  a  programme  of  measurements  was  generated  with  an  emphasis  on 
low-frequency  reverberation;  in  addition,  as  the  activated  towed  array  concept  is 
new  the  desirability  of  collecting  operational  data  under  typical  environmental  con¬ 
ditions  is  urgently  required. 

The  objective  has  thus  been  to  generate  a  database  of  reverberation  returns  which 
not  only  yields  first  and  second  order  parameters  such  as  reverberation  types,  scat¬ 
tering  strength  values,  spectral  characteristics  and  statistics,  but  also  provides  in¬ 
puts  for  further  studies  on  detection,  classification  and  tracking  in  a  reverberant 
environment . 

The  programme  has  up  to  now  had  a  strong  experimental  bias.  Considerable  effort 
has  been  devoted  to  collecting  a  sizeable  and  reliable  database  for  reverberation 
returns  for  a  low-frequency  activated  towed  array  under  a  variety  of  operational 
and  environmental  conditions.  The  initial  data  analysis  has  led  to  seafloor  back- 
scattering  being  identified  as  the  dominating  source  in  the  reverberation  process; 
in  addition,  bottom  backscattering  strength  data  for  typical  areas  have  been  eva¬ 
luated.  An  associated  statistical  analysis  indicated  that  reverberation  for  this  class 
of  sensors  tends  to  be  locally  covariance  stationary,  and  therefore  logarithmic  nor¬ 
malization  techniques  can  be  usefully  applied. 

The  report  commences  with  a  description  of  the  experimental  set-up  and  of  the 
data  processing.  A  review  of  the  theory  of  reverberation  and  of  the  principles  of 
computer  prediction  models  with  a  reverberation  capability  follows,  and  then  the 
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current  theoretical  approaches,  and  related  formulas  for  estimating  reverberation 
characteristics  are  summarized.  In  the  presentation  of  the  results  typical  experi¬ 
mental  data  are  given  in  the  form  of  real-time  colour  displays  of  recorded  intensity 
time  series  of  beamformed  outputs  for  various  waveforms.  A  preliminary  analysis  is 
made  of  the  relevance  of  bottom  backscattering  and  of  the  statistics  of  the  intensi¬ 
ties  of  both  ambient-noise  dominated  data  and  reverberation-dominated  data.  The 
conclusions  of  the  analysis  are  recapitulated  and  summarized  in  the  final  section. 
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2.  The  experimental  set-up  and  data  processing 

Monostatic  reverberation  measurements  were  obtained  by  transmitting  acoustic 
pulses  from  a  high-power  towed  sound  source,  and  receiving  returns  on  the  SA¬ 
CLANTCEN  Prakla-Seismos  horizontally  towed  linear  array.  Figure  1  shows  the 
system  layout  for  reverberation  measurements.  Transmission  and  acquisition  were 
performed  while  SACLANTCEN  R/V  Maria  Paolina  G.  (MPG)  steamed  along 
stable,  straight  courses,  with  a  constant  speed  of,  usually,  5  kn,  along  a  variety  of 
headings. 

The  arrangement  used,  although  not  optimal  for  scattering-parameter  estimation, 
was  selected  as  the  best  way  of  generating  results  close  to  the  reverberation  that 
would  be  measured  under  the  operational  conditions  for  this  type  of  sensor. 


2.1.  GEOGRAPHICAL  AREAS  AND  ENVIRONMENTAL  DATA 

The  trials  were  conducted  in  the  Mediterranean  Sea  at  the  geographical  locations 
indicated  in  Fig.  2.  Two  general  areas  were  involved:  the  Balearic  Basin  and  the 
South  Tyrrhenian  Sea. 

•  The  Balearic  Basin  is  a  region  typical  of  the  Mediterranean  Abyssal  Plains 
and  features  a  rather  homogeneous  seafloor  structure.  A  description  of 
the  geophysical  and  geological  characteristics  of  the  basin  can  be  found  in 
Ref.  [1],  The  top  layers  of  the  seabed  are  made  up  of  clay,  and  clay/silt/sand 
combinations.  Core  samples  previously  collected  in  the  area,  see  Fig.  3  from 
Ref.  [2],  show  also  the  presence  of  top  layers  of  sand.  The  bottom  depth  is 
very  slowly  changing  over  long  distances,  with  rates  of  less  than  1  m  per  mile 
in  the  measurement  region,  and  an  average  water  depth  of  about  2770  m; 
this  makes  the  area  convenient  for  comparison  with  model  predictions.  (The 
presently  available  reverberation  prediction  models  do  not  allow  for  a  varying 
bottom  depth.) 

•  The  South  Tyrrhenian  Sea  is  also  a  deep  water  basin,  with  a  maximum  water 
depth  of  about  3500  m,  but  it  features  a  complex  seafloor  structure  with 
several  peculiarities  such  as  sea  mounts  and  pinnacles,  and  several  islands 
all,  of  volcanic  origin. 

The  measurements  were  conducted  under  winter  and  summer  environmental  con¬ 
ditions  in  order  to  assess  velocity  profile  effects  on  reverberation  returns.  Figure  4 
shows  sound  speed  profiles  over  3000  m  of  depth,  with  an  expanded  view  over  the 
first  500  m.  The  data  relate  to  early  April  and  early  September  and  thus  are  reaso- 
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nabiy  typical  of  average  winter  and  summer  conditions.  Three  shallow-depth  XBTs 
(T4  type,  to  a  depth  of  400  m)  were  made  per  day. 

As  is  known,  in  summer  the  Mediterranean  Sea  features  a  pronounced,  relatively 
shallow  sound  channel  with  a  minimum  around  85  to  150  m.  This  channel  is 
surmounted  by  a  mixed  layer  from  the  surface  to  a  depth  of  about  50  m,  sometimes 
exhibiting  small  surface  ducts.  For  this  situation,  surface  convergence  zones  for 
a  source  and  receiver  in  the  axis  of  the  sound  channel  are  approximately  spaced 
35  km  apart,  about  half  the  spacing  for  an  Atlantic  summer  environment.  In  winter 
the  conditions  are  characterized  by  a  negative  velocity  gradient  which  gives  rise  to 
an  upward  refracting  profile,  although  a  minor  surface  layer  may  be  present. 


Measurements  were  made  with  a  pretty  well  distributed  set  of  wind  speeds  and 
sea-states,  ranging  up  to  sea-state  6  (with  wind  gusts  of  up  to  45  kn).  Sea-state  5 
or  6  is  the  limit  for  safe  operations,  for  the  launching  and  retrieval  of  the  projector. 


2.2.  SEA  TRIAL  CHRONOLOGY  AND  TRANSMITTER  DATA 

Table  1  gives  a  summary  of  the  chronology  of  the  Sea  Trials.  The  table  also  provides 
information  on  the  features  of  the  various  sources  deployed,  since  the  source  is  the 
most  important  element  in  reverberation  measurements. 


TABLE  1 

Chronology  of  sea  trials  and  transmitter  data1 


Date  June  ’83  April  '84  Sept.  '84  Sept.  ’85 

Duration  (days)  4  4  1  5 

Area1  ST  BB  ST  ST-BB 


Type3 

HX-90R 

f.t. 

f.t. 

f.t. 

Configuration4 

single 

b.b. 

b.b.  | 

r  vertical  dipole 
[1.92  m  spacing 

Level  (dB/jtPa/lm) 

203 

209 

209  | 

r  215  (at  370  Hz) 

[200  (at  740  Hz) 

Pattern 

omni 

omni 

omni  | 

f dipole  (at  370  Hz) 
[omni  (at  740  Hz) 

Frequency  (Hz) 

333 

370 

370 

370-740 

1  Depth:  100  m  (nominal). 

1  ST:  South  Tyrrhenian;  BB:  Balearic  Basin. 
3  f  t.:  fiextensional  transducer. 

1  b.b.:  back-to-back. 
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A  detailed  description  of  the  history  of  the  evolution  of  the  transmitters  can  be 
found  in  Ref.  [3j.  Initially  a  single  HX-90R  that  resonated  at  333  Hz  was  available, 
and  this  was  followed  by  a  source  made  out  of  two  flextensional  transducers  connec¬ 
ted  in  a  back-to-back  arrangement  and  resonating  at  440  Hz.  This  configuration, 
although  providing  a  higher  source  level,  was  nevertheless  solely  omnidirectional. 
During  the  September  ’85  Sea  Trial,  a  new  towed-body  sound  source  was  intro 
duced,  still  with  two  flextensional  transducers,  but  along  a  vertical  axis  with  a 
separation  of  1.92  m.  In  the  elevation  plane  the  transmission  pattern  of  this  source 
corresponds  to  a  dipole  with  approximately  -20  dB  sidelobes  at  a  measurement 
frequency  band  centred  around  370  Hz.  This  same  configuration  was  adso  used  for 
measurements  in  a  high  frequency  band  (740  Hz  centre  frequency);  in  this  case  the 
pattern  was  essentially  omnidirectional  except  for  a  null  at  ±  35°  elevation.  The 
cable  length  was  such  that  the  source  trailed  the  towship  by  about  450  m. 


2.3.  DATA  ACQUISITION  SYSTEM  CONFIGURATION 

Monostatic  reverberation  returns  were  acquired  on  the  Prakla-Seismos  towed  array. 
The  array  was  kept  at  a  nominal  depth  of  100  m  and  at  a  horizontal  separation 
from  the  towship  of  900  m.  The  hydrophone  configuration  consisted  of  a  32-element 
array  at  both  frequencies  (370  and  740  Hz);  however  the  hydrophone  spacing  was 
1.96  m  (half-wavelength  frequency  at  382  Hz)  for  measurements  at  370  Hz,  and 
0.98  m  (half- wavelength  at  765  Hz)  for  measurements  at  740  Hz.  Sensors  monitor 
the  array  head  and  tail  depth  and  orientation.  A  simplified  block  diagram  of  the 
data  acquisition  system  is  given  in  Fig.  5.  The  hydrophones’  preamplified  output 
voltages  are  relayed  in  parallel  to  the  acquisition  and  processing  system  aboard 
the  towship,  where  the  signals  are  amplified,  filtered,  digitized  and  beamformed; 
they  are  also  recorded  in  digital  format,  on  a  Bell  and  Howell  High  Density  Data 
Recording  (HDDR)  unit  for  subsequent  off-line  analysis. 

The  Signal  Conditioning  Unit  (SCU)  performs  analog  filtering  and  amplification  (in 
discrete  6-dB  steps),  in  parallel,  on  the  hydrophone  voltages.  The  filter  bandwidth 
is  selected  to  prevent  aliasing  in  the  following  analog-to-digital  (A/D)  conversion. 
The  SCU  gain  was  adjusted  with  a  trade-off  between  maximum  sensitivity  and 
linearity  before  filtering.  A  bank  of  programmable  passband  filter-amplifiers  was 
added  from  the  September  ’85  Sea  Trial,  which  allows  an  increase  in  gain  within 
the  (narrow)  signal  bandwidths. 

The  analog  voltages  are  digitized  by  a  bank  of  11-bit  A/D  converters  in  the  Front 
End  Unit.  The  sampling  frequency  is  chosen  to  be  at  least  three  times  the  maximum 
transmitted  frequency.  The  digital  samples  are  stored  on  the  HDDR  unit  for  off¬ 
line  processing.  For  real-time  processing,  the  data  are  fed  to  WARP  (i.e.  the  Wide 
Application  Real-Time  Processing)  system. 


The  WARP  consists  of  an  interpolator,  a  beamformer  and  a  CSPI  MAP  (Macro 
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Array  Processor)  with  a  fast  one-megaword  RAM,  called  the  DataRam.  These 
devices  are  supervised  by  a  Hewlett-Packard  HP-1000  minicomputer  system  with 
standard  disk  and  magnetic  tape  peripherals.  This  chain  performs  time  domain 
beamforining  and  matched  filtering  on  the  hydrophone  time  series. 

The  interpolator  consists  of  two  cascaded  banks  of  parallel  linear-phase  program¬ 
mable  filters.  The  first  stage  is  actually  used  to  band-pass  filter  the  input  data.  A 
typical  filter  response,  including  finite  word  length  quantization  errors,  is  shown  in 
Fig.  6a.  The  second  stage  does  up-sampling  for  beamforining,  by  zero  padding  and 
low-pass  filtering. 

No  direct  complex  demodulation  is  applied  to  the  input  data;  instead,  a  sequence 
of  filtering,  up-sampling  and  decimating  operations  is  used  to  translate  the  signal 
down  to  baseband.  This  requires  the  various  sample  rates  to  be  congruent  (i.e.  all 
integer  multiples  of  one  another).  The  driving  parameter  is  the  bandwidth  of 
interest,  which  presently  has  a  value  of  35  Hz  at  370  Hz  and  40  Hz  at  740  Hz 
(as  a  consequence  mainly  of  projector  limitations).  As  a  result,  the  final  sample 
rates  are  respectively  70  and  80  Hz. 

To  accurately  form  beams  with  a  fine-grain  angular  resolution,  a  high  sample  rate 
must  be  selected  at  the  beamformer  input.  Typically,  a  resolution  of  0.1°  in  azimuth 
requires,  at  370  Hz,  a  sample  rate  of  about  13  kHz.  The  interpolator  is  therefore 
used  to  raise  the  input  frequency,  which  is  relatively  low  due  the  constraints  imposed 
by  the  A/D  converters  and  HDDR.  In  addition,  it  is  also  used  to  prefilter  data 
in  order  to  allow  decimation  at  the  beamformer  output  as  a  mean  of  avoiding 
aliasing.  Consequently,  the  data  rate  at  the  beamformer  output  is  substantially 
reduced.  A  further  decimation  of  the  various  beam  time-series  takes  place  in  the 
array  processor,  after  which  matched  filtering  and  spectral  analysis  are  applied. 

The  schematic  below  shows  the  sample  rates  at  the  various  stages,  for  both  fre¬ 
quency  bands;  Fig.  6  b,c  shows  the  shape  of  the  anti-aliasing  filters  in  the  interpo¬ 
lator  and  in  the  array  processor,  for  the  low-frequency  band  (around  370  Hz).  This 
selection,  which  results  in  a  set  of  beam  time  series  sampled  at  the  Nyquist  rate, 
is  determined  by  a  number  of  trade-offs  related  to  various  practical  factors  such  as 
device  dynamic  range  and  performance  limitations  and  word  length  and  number  of 
coefficients. 
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The  beamformer  performs  shift-and-add  time  domain  beamforming  on  the  hy¬ 
drophone  inputs,  which  generates  parallel  beam  output  time  series,  with  optional 
shading.  The  beamforming  arrangement  that  was  most  frequently  used  consists 
of  16  beams  panoramically  pointing  into  the  left/right  ambiguous  180°  horizontal 
space,  from  front  to  rear  endfire.  Hamm  shading  is  applied  to  all  the  beams;  this  is 
to  reduce  sidelobes  at  the  expense  of  widening  the  main  beam  (by  about  one  and  a 
half  times).  At  broadside,  the  3-dB  beamwidth  is  thus  about  5.4°.  A  typical  beam 
configuration  as  determined  by  the  conical  symmetry  of  a  linear  array  is  shown 
at  Fig.  7  (horizontal  view  and  side  view.  The  beam  arrangement  is  such  that  it 
keeps  beam  overlap  as  constant  as  possible.  Two  representative  beam  patterns,  at 
broadside  and  at  endfire,  are  given  at  Fig.  8. 


2.4.  WAVEFORMS  AND  SIGNAL  PROCESSING 

In  order  to  adequately  sample  the  reverberation  scattering  function,  continuous 
wave  (CW)  and  frequency  modulated  (FM)  waveforms  were  employed  to  investigate 
reverberation  characteristics  signals  having  good  range  and  good  doppler  resolution 
capabilities. 

As  a  practical  example,  Table  2  lists  the  characteristics  of  typical  transmitted  si¬ 
gnals.  The  codenames  shown  refer  to  the  real-time  colour  displays  of  Sect.  4. 


TABLE  2 

Parameters  of  typical  transmitted  waveforms 


Codename1  (at  fo)  j 
Type 

Pulse  length 

Swept  bandwidth  (FM  signal) 
Repetition  rate 


Ml  (at  370  Hz) 
M7  (at  740  Hz) 
CW  +  FM 
8  s  /  8  s 
10  Hz 
240  s 


CF  (at  370  Hz) 
F7  (at  740  Hz) 
CW  +  FM 
1  s  /  1  s 
10  Hz 
160  s 


¥ 


1  Ping-to-ping  frequency  diversity  was  applied. 


The  following  remarks  must  be  noted: 

(a)  Use  is  made  of  'composite'  signals,  whereby  each  transmission  consists  of 
two  consecutive  sub-pulses  (at  separate  frequencies):  a  CW  pulse  followed 
by  a  linear  FM  pulse  of  equal  duration.  This  allows  direct  comparison  of  the 
scattering  characteristics  of  the  environment  when  wt"eforms  with  contra¬ 
sting  capabilities  (high  doppler,  poor  range  resolution  for  the  CW  pulse;  vice 
versa  for  the  FM  pulse)  are  simultaneously  propagating  into  the  medium. 
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(b)  Long-duration  signals  are  exploited  because  long-range  detection  is  one  of 
the  primary  aims  of  an  active  sonar  system  based  on  the  activated  towed 
array  concept.  This  makes  ii  more  difficult  to  separate  the  various  scattering 
processes. 

(c)  Low  repetition  rates  and  ping-to-ping  frequency  diversity  are  required  in 
order  to  avoid  multiple-time-around  (or  range)  ambiguities.  The  range  am¬ 
biguities  associated  with  major  bathymetric  features  and  coastlines  are  a 
severe  problem  for  a  high-power  active  sonar,  particularly  if  they  lie  in  the 
sound  channel  in  relatively  small  basins  (up  to  a  few  hundreds  miles  across), 
such  as  those  in  the  Mediterranean.  It  is  because  low  repetition  rates  do 
not  completely  overcome  the  problem  that  use  in  also  made  of  frequency 
diversity,  whereby  each  successive  ping  is  transmitted  on  a  different  center 
frequency,  within  the  available  bandwidth  of  the  transducer.  The  combined 
effect  is  elimination  of  the  range  ambiguities,  as  is  apparent  from  the  real 
time  coloured  displays  discussed  in  Sect.  4. 

On  reception,  the  array  processor  applies  conventional  matched  filtering  to  the 
beamformer  outputs.  Beamformer  data  are  demultiplexed  into  the  constituent 
beam  time  series,  which  are  then  individually  filtered  and  correlated.  The  selection 
of  the  processing  parameters  is  such  as  to  adequately  sample  the  delay /doppler 
ambiguity  plane  and  is  made  according  to  Ref.  [4], 


OW  processing.  CW  processing  evaluates  the  spectral  time  history  of  the  indi¬ 
vidual  beam  data  (on  a  ping-by-ping  basis)  by  computing  time-overlapped  Hanning 
shaded  power  spectra.  The  spectra  are  obtained  by  means  of  Fast  Fourier  Transform 
(FFT)  algorithms  which  are  taken  for  time  windows  that  have  the  same  length  as 
the  transmitted  signal  duration  and  are  zero-padded  to  the  appropriate  length  and 
evaluated  over  data  segments  that  have  50%  overlapping  in  time. 

Towship  motion  compensation,  or  ou’n  doppler  nulling  (ODN),  is  applied  to  the  C'W 
data  by  realigning  the  received  power  spectra  as  a  function  of  the  beam  pointing 
direction  and  the  ship's  velocity  vector,  according  to  the  formula 

t  2V,hip 

/d  =  —J—  cos  7 ,  ( 1 ) 

where  fa  is  the  doppler  frequency  (Hz),  monostatic  case;  V.hip  is  the  ship  speed 
(m/s);  X  is  the  acoustic  wavelength  (m);  7  is  the  angle  between  ship  velocity  vector 
and  beam  pointing  direction. 

By  applying  a  correction  equal  to  -  fa  all  returns  from  non-moving  targets  which  are 
located  along  the  maximum  response  axis  (MRA)  of  each  be,-m  pointing  direction 
will  have  spectra  centered  at  zero  doppler  shift.  This  is  true  regardless  of  the  actual 
location  of  the  scatterer  on  the  three-dimensional  MRA  curve,  which  for  a  linear 
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array  is  a  cone  with  its  axis  along  the  array  axis.  (However  the  realignment  is  in 
practice  seldom  perfect  due  to  the  ODN  algorithm  being  unable  to  compensate  for 
ship-motion  fluctuations  (or  deviations)  from  the  indicated  value  and  transmit-to- 
receive  velocity  variations,  with  the  latter  deficiency  being  the  most  significant.) 

The  motion  of  the  towing  platform  thus  introduces,  for  every  beam,  a  relationship 
between  (a)  the  spatial  displacement  of  a  target  with  respect  to  the  beam’s  MRA 
and  (b)  the  doppler  shifted  frequency  of  its  received  echo.  As  will  be  clear  later 
on,  this  is  of  importance  for  returns  from  major  bathymetric  features  and  doppler 
sensitive  waveforms  (long  CW  pulses).  For  a  given  beam  pointing  direction,  returns 
received  from  an  extended  scatterer  (which  may  be  typical  of  reverberation)  throu¬ 
ghout  the  whole  beam  pattern  will  thus  be  affected  not  only  by  amplitude  shading 
but  also  by  a  frequency  shift  as  a  function  of  the  angle  between  the  MRA  and  the 
received  scattering  area. 

For  multiple  beams  formed  with  a  linear  array,  the  frequency  spread  due  to  main 
beam  off-MRA  returns  is  approximately  independent  of  the  beam  pointing  di¬ 
rection.  This  can  be  shown  with  the  following  argument.  Since  ODN  applies  a 
(-2V/X)  cosy  doppler  compensation,  the  resulting  incremental  doppler  shift  of  off- 
MRA,  main-beam  returns  is 

d/d  =  t[cos(7  +  dy)  -  cosy]  ,  (2) 

where  k  =  2V/X  and  dy  is  the  off-axis  angle. 


For  small  dy,  it  is 


d/d  ~  fcdy  siny  ,  (3) 

and  since  for  linear  array  it  is  approximately 

dym„  =  3-dB  beamwidth  ~  y?t,r/  siny  ,  (4) 

where  <fibr  =  3-dB  beamwidth  at  broadside  if  y  is  not  too  close  to  endfire,  and  hence 


d/d  k<fi br ,  (5) 

which  is  independent  of  y.  For  the  experimental  set-up  under  consideration  here, 
Eq.  (5)  gives  d/d  =  0.11  Hz  at  370  Hz  and  d/d  =  0.22  Hz  at  740  Hz. 


LFM  signals  Linear  frequency  modulated  (LFM)  signals  are  processed  through 
a  conventional  quadrature  replica  correlator.  The  de- multiplexed  beamformer  out¬ 
puts  are  correlated  in  two  quadrature  channels  with  a  Hanning  shaded  replica  in 
order  to  reduce  range  sidelobes  effects.  No  doppler  compensetion  is  applied  to  FM 
processing  since  the  maximum  expected  loss  for  the  worst  case,  such  as  stationary 
targets  in  the  endfire  beams,  is  less  than  1  dB  (due  to  the  fact  that  the  modulation 
bandwidths  used  are  small). 
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For  both  the  CW  and  the  LFM  signals,  envelope  detection  is  applied  after  matched 
Altering,  by  squaring  and  adding  the  complex  or  quadrature  channels.  The  data 
are  further  multiplied  by  a  calibration  factor  which  converts  values  to  equivalent 
dBre/iPa  pressure  level  in  the  processed  band,  as  seen  by  a  single  hydrophone. 
This  forms  the  database  for  display  and  further  analysis.  Therefore  the  acquisition 
and  processing  programs  create  large  arrays  of  the  type 

icw  =  /(/., fi./Ji.Ci).  for  CW  data, 

(6) 

$FM  =  for  FM  data, 

where  *hcw  is  the  magnitude  or  complex  valued  processed  CW  output;  is  the 
magnitude  valued  processed  FM  output;  /,  is  the  doppler  frequency  index;  t,  is  the 
time-to-transmission  index;  /3,  is  the  beam  index;  (,  is  the  transmission  number 
(ping)  index. 


To  keep  the  storage  requirements  and  computational  speed  of  these  large  databases 
manageable,  and  yet  not  lose  information,  the  samplings  of  the  time  and  frequency 
indexes  for  the  processed  output  are  selected  such  that  for  the  time  increment  df 
(=f.+i-f.) 


dt  = 


1 

FM  bandwidth 


for  FM  signals. 


(7) 


which  gives  ~  1.5  samples  per  actual  3-dB  resolution  (due  to  the  Hanning  weigh¬ 
ting),  and 

df  =  ^fp  for  CW  signals.  (8) 


where  fp  is  the  CW  pulse  duration  required  by  the  50%  overlap  processing;  the 
resulting  theoretical  overlap  correlation  is  ~  16%  (see  Harris  [5]). 


For  CW  waveforms,  the  frequency  increment  df  (=  /,+ 1  -  /,),  which  is  the  fre¬ 
quency  definition  of  the  spectra,  is  selected  such  that 


in  which 


qit  ? 

2N 


(9) 


where  /,  is  the  processing  sample  rate;  and  N  is  the  nearest  larger  integer  to  obtain 
qi  >  1;  and  2N  is  the  FFT  length.  This  selection  makes  the  frequency  definition 
(or  bin)  ~  2  times  smaller  than  the  theoretical  frequency  resolution  of  the  CW 
waveform  (with  Hanning  shading). 


Having  defined  the  selection  for  the  increments  of  the  various  indexes,  the  criteria  to 
specify  the  sire  of  the  databases  are  based  on  considerations  related  to  the  physical 


-  10  - 


SACLANTCEN  SR- 112 


limitations  for  maximum  transfer  rates.  These  limitations  are  not  described  here, 
but  in  general  the  coverage  in  frequency  is  such  as  to  have  a  fixed  frequency  range 
at  a  given  carrier  irrespective  of  other  waveform  parameters  (for  instance  it  is  4  Hz 
at  370  Hz  centre  frequency). 


» 
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3.  Review  of  reverberation  theory 


3.1.  INTRODUCTION 


Reverberation  is  the  process  that  describes  the  time  variation  of  the  total  scatte¬ 
red  sound  field  observed  at  a  point  of  reception  following  transmission  of  a  signal 
(Ol’shevskii  [6]).  Thus,  it  is  associated  with  the  intrinsic  statistical  inhomogeneity 
of  the  underwater  medium.  Reverberation  is  fundamentally  a  non-stationary  sto¬ 
chastic  process. 

Reverberation  is  commonly  classified  into  two  categories: 

•  Volume  reverberation  due  to  the  scattering  from  inhomogeneities  in  the  ocean 
medium  itself,  such  as  those  attributable  to  biological  origin  (e  g.  the  mi¬ 
croorganisms  of  the  Deep  Scattering  Layer,  for  instance,  and  fish)  or  thermal 
irregularities  or,  possihlv,  internal  waves. 

•  Boundary  reverberation  arising  from  the  scattering  associated  with  the  inte¬ 
raction  of  sound  with  the  medium  discontinuities  represented  by  the  ocean 
boundaries.  This  is  further  classified  into: 

-  surface  reverbi  ration  caused  by  the  scattering  of  sound  at  or  near  the 
sea  surface,  due  to  ocean  surface  waves  or  to  entrapped  air  bubbles, 
and 

bottom  reverberation  due  to  the  scattering  of  sound  by  the  inhomo- 
geneities  of  the  ocean  bottom  and  the  irregularities  of  the  bottom 
surface. 

Two  theoretical  approches  have  been  taken  to  describe  the  statistics  of  the  backscat- 
tered  sound  fields.  Both  approaches,  by  analogy  with  the  Huygens- Fresnel  principle 
of  physical  optics,  consider  the  scattered  sound  field  as  consisting  of  elementary 
waves  in  mutual  phase  interference.  In  the  first  classical  or  'physical'  approach  a 
solution  of  the  appropriate  dynamical  equation  of  the  propagation  is  searched,  in¬ 
troducing  the  stochastic  mechanism  of  random  scattering  as  a  spatial  perturbation 
This  method  incorporates  the  (statistical)  parameters  of  the  scattering  process  into 
the  Helmoltz-Kirchhoff  equations  for  a  lossless  inhomogeneous  medium:  it  therefore 
relates  the  statistics  of  the  scattering  mechanism  with  the  geophysical  bulk  para 
meters  of  the  medium.  However  its  complexity  limits  its  usefulness  in  handling 
complex  geometries,  waveforms  and  higher-order  statistics. 

In  the  second  approach,  the  quasi-phenomenologtcal  model  of  Faure  (7j,  Ol’shev- 
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skii  [6],  and  Middleton  [8j,  a  random  distribution  of  independent  point  scatterers 
in  the  medium  is  introduced  in  the  dynamical  equations  of  propagation.  The  me¬ 
dium  itself  is  treated  as  infinite  and  everywhre  homogeneous  and  isotropic.  First 
order  scattering  is  commonly  assumed:  that  is,  the  scatterers  do  not  mutually 
interact  and  simply  reradiate  part  of  the  incident  field  in  all  directions.  This  is 
statistically  equivalent  to  the  assumption  that  the  scattering  process  is  a  Poisson 
process  independent  both  in  space  and  time.  Ray  theory  is  usually  invoked  to  de¬ 
termine  geometric  boundaries  of  the  illuminated  scatterers  and  the  received  fields; 
additionally,  the  boundary  conditions  can  be  expressed  by  means  of  a  dynamical 
impulse  response  function  which  represents  the  relationship  between  the  incident 
and  reradiated  field.  Although  this  approach  does  not  directly  relate  scattering 
with  physical  medium  parameters,  it  handles  complex  situations  and  higher-order 
statistics. 


3.2.  SCATTERING  STRENGTH  AND  PREDICTION  MODELS 


In  general,  the  mean  reverberation  intensity  (a  function  of  time)  depends  on  the 
distribution  and  the  characteristics  of  the  scatterers  in  the  illuminated  region.  The¬ 
refore  the  time-varying  average  is  related  not  only  to  geometry  and  system-related 
factors  (such  as  source  power  and  effective  pulse  length,  transmit  and  receive  beam 
patterns  and  sound  speed  profiles)  but  also  to  factors  related  to  the  scattering 
mechanism  itself,  such  as  the  mean  number  and  the  acoustical  properties  of  the 
scatterers  within  the  illuminated  area. 

The  latter  factors  are  usually  expressed  by  a  single  parameter,  the  scattering  strength, 
which  is  a  convenient  way  of  quantitatively  describing  the  effect  of  reverberation. 
It  represents  the  ratio  between  the  intensity  of  the  sound  scattered  by  a  unit  area 
or  volume  referred  to  a  distance  of  1  m  and  the  incident  plane- wave  intensity.  This 
ratio  is  commonly  expressed  in  decibels,  and  for  volume  scattering  is 

S„  =101ogs„  .lologf^^)  (10) 

and  for  boundary  scattering 

5„,  =  101ogsb,  ^Ologf^f^).  (11) 


The  scattering  strength  is  thus  defined  in  a  similar  manner  to  the  target  strength.  It 
is  used  in  computations  of  the  echo  to  reverberation  level,  and  three  types  of  scatte¬ 
ring  strengths  can  be  distinguished  according  to  the  different  types  of  reverberation 
(5,  for  volume,  5,  for  surface  and  Sb  for  bottom  scattering). 
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By  assuming  the  scattering  strength  constant  over  the  incremental  area  there  is  con¬ 
siderable  simplification  of  the  propagation  equations;  nevertheless,  most  situations 
of  practical  interest  are  still  too  complex  to  permit  a  closed-form  solution  to  the 
problem  of  estimating  the  received  intensity  and  thus  it  is  necessary  to  apply  nu¬ 
merical  integration  techniques.  Consequently,  one  resorts  to  computer  models  that 
include  a  reverberation  prediction  capability,  such  as  the  GENERIC  Sonar  Model 
(Weinberg  [9]),  NISSM  (Weinberg  (12|),  the  NRL  Reverberation  Model  (Pranchi 
et  al.  (11])  and  LIRA  (Hoffmann  [10]). 

The  method  used  by  the  computer  models  to  estimate  reverberation  levels  is  rather 
straightforward  and  relies  on  linear  and  first-order  approximations  of  the  scat  tering 
process;  that  is,  no  allowance  is  made  for  secondary  scattering  effects  (scattering 
from  scattering)  or  mutual  coupling,  and  the  scattering  strength  over  the  illumina¬ 
ted  incremental  area  is  constant . 

The  received  level  at  a  sensor  following  transmission  of  a  signal  of  power  P  and 
effective  duration  fp  is  evaluated  by  considering  the  closed  ray  path  that  is  formed 
between  the  source  and  an  incremental  scatterer,  and  another  to  the  receiver.  The 
acoustic  pressure  at  the  end  of  this  pair  of  eigenrays  is  given  by 

/„,  =  P.D.rxSrdtrfoD,  ,  (12) 

where  P,  is  the  source  power;  l,„  is  the  received  intensity;  D.  is  the  source  direc¬ 
tional  gain  in  the  direction  of  the  outgoing  ray,  D,  is  the  receiver  directional  gain 
in  the  direction  of  the  incoming  ray;  r/t  is  the  propagation  loss  in  the  forward  path; 
ijr  is  the  propagation  loss  in  the  backward  path;  sr  is  the  scattering  strength;  d^ 
is  the  incremental  scattering  area  or  volume. 

Since  the  effective  duration  of  the  signal  is  fp,  the  contribution  to  the  received  level 
comes  front  all  forward  and  received  (multipath)  eigenrays  which  satisfy 

to  tp<ff  +  f,  <t0,  (13) 

where  to  is  the  reference  time  at  which  the  power  is  measured  and  1 1  and  t,  are  the 
forward-  and  receive-path  travel  times  respectively. 

Since  the  resultant  mean  reverberation  pressure  level  is  the  random-pnase  addition 
ol  the  individual  contributors,  it  is 

P}„  -  f  f°  |/r„!:  df  <\<t> ,  (14) 

J*  Jt o  If 

where  ♦  is  the  illuminated  area  or  volume. 


In  the  evaluation  of  the  received  level  account  must  be  taken  of  ail  possible  paths 
from  the  source  to  the  scattering  region  and  back  to  the  receiver.  The  evaluation 
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of  Eq.  (14)  is  accomplished  by  a  numerical  summation  over  a  ‘reasonable'  number 
of  significant  paths  and  over  the  illuminated  area.  Equation  (14)  also  shows  that 
it  may  be  difficult  to  separate  contributions  coming  from  different  spatial  regions 
(which  are  thus  of  different  types)  if  the  travel  times  of  their  eigenray  pairs  fall 
within  the  same  window;  this  comment  is  particularly  applicable  to  long-range 
reverberation. 

The  practical  complexity  of  this  solution  shows  how  necessary  a  computer  model 
is  for  the  prediction  of  the  received  levels  and  the  prediction  of  the  scattering 
parameters.  However  a  simplified  approach  that  is  useful  for  first-order  parameter 
assessment  can  use  the  standard  sonar  equation  (Uriclt  [13] )  in  the  form 

RL  =  SL  -  TLj  +S,C  +  101ogilt,  (15) 

where  RL  is  the  received  power  level  (dBre/iPa),  SL  is  the  source  level  (dBre/iPa), 
TLj  is  the  two-way  propagation  loss  (dB),  S,c  is  the  scattering  strength  (dB), 
is  the  integral  over  the  scattering  region. 

Such  a  simplification  is  applicable  when  the  various  elements  of  Eq.  (15)  are  se¬ 
parable;  and  it  is  convenient  to  use  it  when  the  factors  can  be  easily  expressed, 
as  can  for  instance  the  propagation  loss,  by  means  of  an  analytical  spreading  law. 
However,  the  equivalent  scattering  region  4„  must  be  estimated  as  a  function  of 
the  transmission  parameters  and  of  the  source  and  receiver  beam  patterns. 


3.3.  VOLUME  vs  BOUNDARY  REVERBERATION  AT  LOW  FREQUENCY 


Volume  reverberation  is  mostly  of  biological  origin,  and  is  commonly  associated 
with  the  Deep  Scattering  Layer;  therefore  volume  scattering,  at  the  low  end  of  the 
frequency  spectrum  (down  to  2  kHz),  is  most  likely  caused  by  resonance  effects  from 
the  gas-filled  swimbladders  of  the  larger  types  of  fish.  At  low  enough  frequencies, 
below  1  kHz,  where  wavelengths  are  in  the  order  of  more  than  1  m,  volume  scattering 
strength  can  therefore  be  expected  to  be  small. 

Measurements  of  S,  in  the  Mediterranean  Sea  (Doutt  [14],  lowest  frequency  analy¬ 
sed  1.6  kHz)  showed  the  spectral  behaviour  as  peaking  at  frequencies  in  the  range 
of  4  to  8  kHz.  As  the  measured  fall-off  rate  of  S,  vs  frequency  at  the  lower  end 
was  close  to  the  theoretical  40  dB/decade  expected  from  the  resonant  nature  of 
the  scattering  from  swimbladders,  this  ruled  out  the  possibility  of  large  scatterers 
(swimbladder  diameters  larger  than  a  few  tens  of  cm)  having  significantly  contri¬ 
buted  to  the  scattering  process.  A  value  of  -110  dB/m  for  the  volume  scattering 
strength  has  been  taken  for  computer  model  predictions.  This  is  a  conservative 
figure  derived  by  extrapolating  from  data  for  higher  frequencies  with  the  scattering 
taking  place  at  the  depth  of  maximum  reverberation. 
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Although  the  volume  scattering  strength  is  much  lower  than  that  of  the  bounda¬ 
ries,  it  may  nevertheless  become  predominant  at  long  ranges  because  the  ensonified 
volume  may  grow  large  enough  with  respect  to  the  ensonified  areas.  The  ratio 
between  volume  and  boundary  reverberation  at  a  given  distance  can  be  approxi¬ 
mately  expressed,  from  Eq.  (15),  assuming  all  other  factors  identical,  as 

RL„  r  RLbr  =  5,  -  Sb  +  (lOlogJ  *(yj,i>)dl-'-  lOlog  J  *(*.!*)  d.d  ).( 1«| 

where  RL»  r  is  the  received  level  due  to  volume  reverberation,  RLb  r  is  the  received 
level  due  to  boundary  reverberation,  5,  is  the  column  scattering  strength.  Si,  is  the 
boundary  scattering  strength,  is  the  beam  pattern  (function  of  azimuth 

and  elevation  angles),  T  is  the  scattering  volume  at  range  R ,  A  is  the  scattering 
area  at  range  R. 

The  usual  simplifying  assumptions  for  the  evaluation  of  the  scattering  area  and 
volume  show  that  for  an  unbounded  medium  (which  is  an  unrealistic  hypothesis 
for  long  ranges  where  the  elevation  angles  are  limited  by  the  water  depth)  the 
magnitude  of  the  expression  within  parenthesis  in  Eq  ( lb)  grows  at  approximately 
the  rate  of  lOlogff.  Now  since  the  ratio  between  volume  and  boundary  scattering 
strength  can  at  these  frequencies  confidently  be  expected  to  be  at  least  -50  dB  or 
larger,  it  foUows  that  volume  reverberation  may  exceed  boundary  reverberation  at 
ranges  of  100  km  or  more.  However,  since  reverberation  decays  as  a  function  of 
range,  volume  reverberation  should  exceed  boundary  reverberation  for  moderately 
high  source  levels  that  are  well  below  the  ambient  noise  levels 

This  is  illustrated  at  Fig.  9  which  shows  the  GENERIC  Sonar  Model  prediction  of 
the  reverberation  levels  for  a  system  configuration  and  environment  simulating  the 
present  experimental  one.  Source  and  receiver  are  located  in  the  sound  channel, 
reverberation  is  estimated  for  a  transmission  of  a  CW  pulse  of  8  s  at  370  Hz  and 
a  source  level  of  213  dBrepPa®  1  m.  Reception  at  the  broadside  beam  of  a  32- 
element,  shaded,  towed  array  is  simulated.  Commonly-used  values  for  the  surface 
and  bottom  scattering  strengths  have  been  selected  from  the  Chapman- Harris  and 
McKenzie  formulas  respectively.  Volume  scattering  data  have  been  taken  equal  to 
-110  dB/m  at  the  usual  deep  scattering  layer  (DSL)  depths,  as  postulated  above. 
The  levels  for  the  three  forms  of  reverberation  show  that  volume  reverberation 
can  be  expected  to  be  some  20  dB  below  the  other  forms  of  reverberation  over 
several  tens  of  kilometers  Volume  returns  may  eventually  become  larger  at  those 
ranges  where  reverberation  levels  are  well  below  noise  levels.  These  are  the  reasons 
why  volume  reverberation  effects  have  been  neglected  in  the  further  analysis  of  the 
present  data. 


J.4.  SURFACE  REVERBERATION  THEORY 

The  theory  of  the  scattering  of  sound  from  periodic  or  random  surfaces  is  well 
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established  (see  Fortuin  (15] )  with  the  basic  work  of  Lord  Rayleigh  [Id]  and  the  ex¬ 
tensions  of,  among  others,  Marsh  [17],  and  the  approach  of  Eckart  [19]  and  Brekhov- 
skikh  [20],  which  are  based  on  the  Kirchhoff  approximation  to  the  Helmholtz  in¬ 
tegral  of  the  scattered  field  induced  at  the  boundary  by  the  incident  wave.  The 
composite-roughness  model  for  the  diffraction  effects  due  to  the  surface  is  an  ap¬ 
proach  which  presently  receives  widespread  acceptance  and  essentially  merges  re¬ 
sults  from  both  the  Rayleigh-based  methods  with  those  based  on  the  Kirchhoff 
approximation  (see  McDaniel  [21]).  Hence,  a  geometrical  optics  approximation  is 
applicable  (Eckart)  if 

a,  =  Jthsinfl  >  1  ,  (17) 

where  a,  is  the  Rayleigh  roughness  parameter;  k  =  2r/A  is  the  acoustic  wavenum¬ 
ber;  A  is  the  acoustic  wavelength;  h  is  the  rms  surface  height  (surface  roughness 
parameter);  9  is  the  grazing  angle  at  the  surface.  And  this  condition  (17)  applies  for 
the  high-frequency  situation  as  well  as  for  forward  near-specular  or  monostatic  high 
grazing  backscatter.  This  approximation  predicts  that  the  scattered  field  angular 
intensity  is  dependent  on  the  local  specular  reflection  associated  with  the  surface 
slope  normal  to  the  impinging  sound.  As  a  consequence  the  scattering  process  is,  to 
a  large  degree,  a  coherent  phenomenon,  with  the  intensity  dependent  on  the  mean 
square  slope  of  the  surface  but  independent  of  wavelength. 

If  however  ot  <  1  as  would  be  the  case  for  low  acoustic  frequencies  impinging  upon 
the  surface  at  small  grazing  angles,  then  linear  scattering  theory  (Rayleigh-Marsh) 
should  provide  a  good  qualitative  description  of  the  scattered  field  from  a  complex 
boundary  ensonified  by  a  plane  acoustic  wave.  This  description  considers  the  surface 
acting  like  a  diffraction  grating,  and  therefore  the  scattering  is  essentially  resonant. 
For  a  purely  periodic,  infinite  surface,  scattering  takes  place  in  discrete  directions, 
and  at  discrete  frequencies  given  by 

sind,  -  sin  d,  +  nA/A  ,  (18) 

where  t),  and  9-,  are  respectively  the  scattered  and  incident  angles  measured  from 
the  vertical  direction  (see  Fig.  10);  A  and  A  are  respectively  the  acoustic  and  the 
surface  wavelengths;  and  n  is  the  scattering  order.  For  n  =  0  scattering  is  the 
‘specular’  reflection,  and  is  a  coherent  process  at  the  frequency  of  the  incident 
wave;  n  =  1  gives  the  first-order,  diffuse  scattering  component  which  is  shifted  in 
frequency.  The  theory  shows  that  the  zero  order  and  the  first  order  are  the  most 
relevant  contributors  to  the  amplitude  of  the  scattered  field.  Equation  (19)  reflects 
the  resonant  nature  of  the  scattering  process,  since  a  given  acoustic  wavelength 
will  ‘select’  a  particular  wavelength  of  the  spectrum  of  the  surface  roughness  (or  a 
narrow  band  thereof).  The  diffuse  re-radiation  will,  however,  not  be  truly  discrete 
due  to  the  randomness  and  finite  size  of  the  illuminated  surface.  Marsh,  Shulkin 
and  Kneale  [18])  have  extended  the  analysis  to  random  surfaces,  and  have  computed 
the  values  of  the  first-order  scattered  field,  which  turns  out  ti  be  a  function  of  the 
surface  wave  spectrum  only.  The  amplitude  of  the  first-order  backscattered  field  is 


-  17 


SACLANTCEN  SR- 112 


approximately  given  by 

At  -  4*in40PV.)  (19) 

ff* 

where  B  is  a  dimensionless  constant  and  Pz(u/„ )  is  the  power  spectrum  of  the  sea 
surface  elevation  vs  wave  frequency.  If  P2(ww)  is  assumed  equal  to  the  Pierson- 
Moskowitz  sen-state  spectrum  assumption  (see  also  Fig.  11)  then 

PV.)  =  r92wwsexp(-0.7V/u>tt;4),  (20) 

where  g  is  the  acceleration  of  gravity,  F  is  a  normalization  factor,  and  v  is  the 
wind  speed.  By  neglecting  the  exponential  term  in  Eq.  (20)  (for  a  fully  developed 
sea),  ill  from  Eq.  (19)  should  be  independent  of  and  be  dependent  only  on 
geometrical  factors. 


These  conclusions  have  several  consequences.  Firstly,  the  amplitude  of  the  scatte¬ 
red  field  should  be  reasonably  independent  of  the  sea- state  if  the  acoustic  frequency 
involved  is  above  roughly  SO  Hz,  because  the  selected  wave  frequency  will  be  ap¬ 
proximately  independent  of  sea-state.  However,  there  should  be  a  dependence  on 
the  sea-state  direction  due  to  the  geometrical  diffraction  mechanism  itself.  Se¬ 
condly,  since  the  sea  surface  is  a  time-dependent  moving  boundary,  the  motion  of 
the  ocean  waves  should  impose  a  frequency  shift  on  the  first  order  scattered  field. 
This  doppler  effect  can  be  estimated  by  considering  the  relationship  between  sur¬ 
face  wavelength  A  and  wave  frequency  /..  For  gravity  waves  in  deep  water  the 
surface  wavenumber  is  given  by 


K  = 


(21) 


where  w„  =  2ir By  combining  Eqs.  (18)  and  (21)  and  for  the  predominant  first- 
order  scattering  (n  =  1)  and  assuming  i),  ~  -i?j  close  to  grazing  incidence  (d  ~ 
90°),  which  is  the  case  for  reverberation  from  ranges  where  R,  3>  d, ,  (in  which  d, , 
is  the  source-receiver  depth  and  R,  is  the  range  from  which  surface  reverberation 
is  received),  then 


where  ft  is  the  doppler  shift  of  the  backscattered  field  and  appears  to  be  related 
only  to  the  incident  wavelength.  This  constant  shift  should  be  ~  0.87  Hz  at  370  Hz 
(acoustic  frequency)  and  ~  1.24  Hz  at  740  Hz.  The  polarity  of  the  doppler  shift  in 
Eq.  (22)  is  determined  by  the  sign  of  the  cosine  of  the  angle  between  the  direction 
of  the  incident  energy  and  the  predominant  sea  direction. 

The  overall  behaviour  of  the  scattering  strength  as  a  function  of  the  grazing  angle  is 
therefore  expressed  by  the  composition  of  two  curves  representing  the  two  regions 
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of  near  specular  and  diffuse  backscatter.  There  should  be  little  or  no  frequency 
dependency  for  high  grazing  angles,  but  it  can  be  expected  in  the  medium  to  low 
grazing  region.  As  far  as  experimental  data  for  low-frequency  surface  scattering 
are  concerned,  the  results  of  Chapman  and  Harris  [22]  for  diffuse  scattering  and 
Chapman  and  Scott  [23]  for  high  grazing  angles  sue  in  widespread  use.  (These  fire 
the  ‘default’  values  of  most  computer  prediction  models.)  For  diffuse  scattering 
(medium  to  low  grazing  angles)  Chapman  and  Harris  fitted  the  following  formula 
to  their  experimental  data  (to  within  2  dB,  on  average,  from  the  measured  values): 

S,  =  3.3«log(0/3O°)  -  42.4  log  k  +  2.6, 

k=158(v/1/3)  ,  (23) 

where  S ,  is  in  dB;  0  is  the  grazing  angle  in  degrees;  v  is  the  wind  speed  in  knots; 
and  /  is  the  acoustic  frequency  in  Hz.  Note  that  Eq.  (23)  predicts  a  sea-state 
dependence  which  1.  not  readily  derived  from  the  Rayleigh- Marsh  theory.  For  high 
grazing  angles,  Chapman  and  Scott  fitted  their  data  to  the  Eckart’s  theory  to  obtain 

S,  -  -101og(8ir£:)  -  2.17c-2/  tan2  8 , 

eJ  =  0.003  +  0.003o.  (24) 


Figure  12a  shows  the  behaviour  of  the  composite  surface  scattering  strength  as 
predicted  by  the  formulas  for  a  number  of  frequencies. 

3.5.  BOTTOM  REVERBERATION  THEORY 

The  theory  of  the  scattering  of  sound  from  the  surface  can  in  principle  be  equally 
well  applied  to  the  scattering  of  sound  from  the  bottom.  However  while  the  surface 
can  be  considered  a  perfect  pressure-release  surface,  it  is  only  approximately,  and 
only  for  very  hard  bottoms,  that  the  seafloor  can  be  considered  a  perfectly  rigid 
boundary.  Seafloor  scattering  is  thus  a  more  complex  process  and  usually  requires 
taking  into  account  the  energy  that  propagates  into  the  sub-bottom  layers,  and  is 
accordingly  much  more  difficult  to  predict,  particularly  at  low  frequencies,  where 
bottom  penetration  is  more  pronounced.  Bottom  reverberation  is  therefore  com¬ 
monly  expressed  as  a  function  of  the  seafloor  geophysical  properties  (sand,  rock, 
clay  or  silt)  or  type,  in  addition  to  the  seabed  roughness. 

A  composite  model  can  also  generally  be  assumed  for  the  scattering  from  the  sea¬ 
floor.  At  medium  and  small  grazing  angles  a  diffuse,  incoherent  form  of  scattering 
as  indicated  by  the  linear  scattering  theory  is  appropriate,  whereas  at  high  grazing 
angles,  a  facet  reflection  model  is  better.  Although  first-order  bottom  scattering  is 
primarily  related  to  the  roughness  of  the  seafloor,  it  is  usually  not  easy  to  determine 
the  bottom  roughness.  A  common  way  of  relating  scattering  to  grazing  angle,  at 


Saclantoen  SR- 1 1 2 


a  given  frequency,  and  in  the  absence  of  a  description  of  the  bottom  structure,  is 
to  assume  that  scattering  is  in  accordance  with  a  Lambert’s  rule  of  physical  op¬ 
tics  for  rough  surfaces,  which  assumes  the  amplitude  of  the  scattered  sound  to  be 
proport  ional  to  the  sine  of  the  grazing  angle.  That  is 

/,  -  p/j  sin  6\  sin  0,  ,  (25) 

where  I,  is  the  intensity  of  scattered  sound;  /,  is  the  intensity  of  incident  sound;  0, 
is  the  grazing  angle  at  incidence;  0,  is  the  grazing  angle  in  the  scattering  direction; 
p  is  the  proportionality  constant. 

Thus  the  scattering  strength  (in  dB)  for  the  case  of  monostatic  backscattering, 
where  0;  ~  6,  is 


Sb(0)  =  10  log 


(!) 


10 log/i  -f  lOlogsin2 


(26) 


Literature  data  in  general  show  a  qualitatively  good  agreement  with  Lambert’s 
rule,  see  for  instance  Schmidt  (24]  and  Urick  (13).  The  behaviour  of  the  bottom 
scattering  strength  vs  grazing  angle  resulting  from  the  composite  model  of  diffuse 
scattering  expressed  by  Lambert’s  rule,  with  p  equal  to  -2?  dB,  and  the  facet 
reflection  model,  is  shown  in  Fig.  12b.  Typically,  experimental  data  show  an  initial, 
rapid  decay  for  near  normal  incidence,  followed  by  a  plateau  region  in  which  Sb 
does  not  strongly  depend  on  the  grazing  angle.  At  very  small  grazing  angles  5b 
should  decay  quickly;  however  measurements  of  the  scattering  strength  at  very  small 
grazing  angles  are  particularly  difficult.  Although  Eq.  (25)  is  rather  empirical,  it  is 
nevertheless  convenient  because  it  requires  the  estimation  of  the  single  parameter 
ft  which  in  general  can  be  taken  to  characterize  a  bottom  type. 

Reported  values  of  low-frequency  bottom  backscattering  strength  indicate  that  it 
is  the  type  of  bottom  more  than  its  roughness  that  is  indeed  the  dominating  factor. 
The  same  data  usually  show  a  small  or  no  frequency  dependence  in  the  sub-kilohertz 
region.  It  is  also  important  to  point  out  that  bottom  backscattering  should  exhibit 
no  frequency  shift,  unlike  surface  scattering,  since  the  seafloor  is  stationary. 
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4.  Experimental  results 

Examples  of  data  sets  of  reverberation  measurements  obtained  with  an  activated 
array  towed  in  deep  water  are  given  in  Figs.  13  to  19.  The  data-sets  are  presen¬ 
ted  in  the  form  of  colour-coded  real-time  displays  of  the  returns  received  following 
the  transmission  of  various  composite  waveforms.  The  system  configuration  is  the 
standard  one  for  a  (quasi  )  monostatic  towed  source/array.  Table  3  provides  a  sum¬ 
mary  of  information  relevant  to  the  runs  which  generated  the  data  in  the  displays. 
Transmission  characteristics  are  deducible  by  referring  to  Table  3,  which  gives  the 
time  of  the  runs  and  also  geographical  and  environmental  information.  The  last 
column  in  the  table  is  a  cross-reference  to  Figs.  20  and  21,  which  schematically 
show  the  positions  of  the  various  runs  with  respect  to  the  basin  boundaries. 


TABLE  3 

Refetence  data  for  the  real-time  colour  displays 


Figure 

no. 

Signal1 

Date/ 
time  (r) 

Area1 

Position 

(mid-run) 

Ship 

heading 

Sea- 

state 

Ref. 

figure3 

13 

C'F 

10  Apr.  ’84, 
11:49 

BB 

41*15' 

07*24' 

N, 

E 

010° 

2 

20- A 

14 

Ml 

10  Apr.  ’84, 
15:27 

BB 

41*35' 

07*31' 

N, 

E 

010° 

2 

20- B 

15 

Ml 

24  Sept.  ’84, 

1  <  :02 

ST 

39*26' 

13*08’ 

N, 

E 

180” 

3 

21- A 

16 

Ml 

4  Sept.  ’85, 
09:37 

ST 

40°  14' 
12°36' 

N, 

E 

o 

© 

r*s 

2 

21-B 

17 

M7 

4  Sept.  ’85, 
10:41 

ST 

40°  16' 
12°29' 

N, 

E 

300° 

2 

21-C 

18 

M7 

6  Sept.  ’85, 
08:06 

BB 

41°48' 
07°  40' 

N, 

E 

010° 

4 

20-C 

19 

Ml 

6  Sept.  '85, 
10:58 

BB 

4  2°  00' 
07°45' 

N, 

E 

010° 

4 

20- D 

1  Refer  to  Table  2. 

3  BB:  Balearic  Basin;  ST:  South  Tyrrhenian 
3  See  text. 


Figure  13  displays  only  the  processed  FM  data  of  the  composite  signal  CF.  The 
figure  shows  a  multi-ping,  multi-beam  real-time  colour-coded  display.  Two-way 
slanl  sonar  range ,  expressed  in  km,  is  displayed  vertically  and  the  various  beam 
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pointing  directions  (refer  for  instance,  for  the  first  two  figures,  to  the  heamformer 
set-up  at  Fig.  7)  are  given  horizontally;  0°  is  forward  endfire,  90°  is  broadside  and 
ISO'  is  rear  endfire  Each  beam/raiige  set  is  made  up  of  sequential-pings  data 
The  received  power  level  is  expressed  in  calibrated  dB  te //Pa  at  the  output  of  the 
matched  filter  and  is  colour-coded  according  to  the  scale  at  the  lower  right  hand  side 
in  the  figure.  Since  echoes  are  recorded  with  a  moving  platform,  non-moving  targets 
generate  ‘tracks’  on  this  multi-ping  display,  with  geometries  that  are  related  to  the 
relative  viewing  angle  from  the  towship.  For  instance,  over  short  time  intervals, 
when  a  target  remains  within  a  single  beam,  a  fixed  target  will  appear  moving 
inward  if  at  the  forward  endfire,  outward  if  at  the  rear  endfire,  and  at  constant 
range  if  at  the  broadside  beam.  Over  longer  time  periods,  a  fixed  target  will  also 
move  from  beam  to  adjacent  beam. 

The  other  figures,  which  pertain  to  the  composite  signals  Ml  (at  370  Hz)  and  M7 
(at  740  Hz),  include  both  FM  and  CW  displays.  The  FM  display  on  the  left  side  is 
identical  to  that  of  Fig.  13.  The  CW  display,  on  the  right  side,  shows  the  received 
power  spectrum  for  each  beatri  as  a  function  of  range,  and  thus  gives  tile  spectra) 
time  history  of  the  CW  returns.  The  various  beam  looking  directions,  which  are 
the  same  as  for  the  FM  display  because  the  returns  are  transmitted,  received  and 
processed  simultaneously,  are  arranged  as  a  cascade  of  smaller  displays  of  received 
intensity  versus  two-way  sonar  range  and  frequency  (doppler  shift).  The  range  (or 
time)  variable,  which  runs  vertically,  covers  the  same  extent  as  in  the  FM  display. 
Ping-to-ping  exponential-decay  averaging  is  applied  to  the  data  and  the  number  of 
pings  used  for  averaging  is  given  by  the  parameter  AV  in  the  figure.  Usually  a  small 
number  of  pings  (3  or  4)  for  the  running  average  is  sufficient  to  smooth  noise  and 
reverberation,  thus  providing  a  reasonably  stable  and  uniform  pattern  that  slowly 
changes  as  a  function  of  the  environment.  As  in  the  FM  case,  the  display  shows 
calibrated  colour-coded  power  levels  at  the  output  of  the  processor.  There  is  no 
need  to  rescale  the  results  since  both  signals  have  the  same  duration  and  hence  the 
same  energy,  and  they  therefore  perform  identically  against  a  white  gaussian-noise 
background.  ODN  is  applied  to  the  CW  data;  all  stationary  scatterers  in  the  main 
beam  of  any  beam  looking  direction  provide  returns  which  are  centered  at  zero 
doppler  shift. 

The  most  prominent  features  of  the  reverberation,  and  its  relationship  with  the 
transmitted  waveform,  are  readily  apparent  from  the  figures,  which  show  the  result 
of  linear  processing  of  the  received  intensit  ies  in  both  space  and  t  ime.  The  ambient 
noise  over  the  period  of  a  run  (usually  about  one  hour)  appears  generally  anisotro¬ 
pic,  and  non-stationary  due  to  towship  noise  contamination  and  to  medium-range 
or  nearby  shipping  contamination.  The  initial  returns  have  very  high  levels,  which 
are  larger  the  closer  the  beam  is  to  broadside,  and  progressively  decay,  depending  to 
some  extent  on  the  waveform  characteristics,  below  the  average  ambient  noise  level. 
LFM  waveforms,  with  higher  range  resolution,  show  that  the  .tructurc  of  the  ini¬ 
tial  returns  is  dominated  by  regularly  spaced  'peaks’,  which  are  associated  with  the 
multiple  bottom-surface  reflections  (the  so-called  fathometer  effect).  Long-range 
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echoes  are  always  easily  related  to  the  backscattering  from  the  basin  boundaries 
(as  a  visual  comparison  of  the  figures  with  the  locations  of  the  runs  confirms)  or  to 
seafloor  highlights  such  as  seamounts.  The  next  section  amplifies  the  data  analysis 
and  gives  further  interpretation  of  the  results. 
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5.  Interpretation  of  results  and  model  comparison 

In  Stibsect.  3.3  it  has  been  shown  that  volume  scattering  did  not  (day  an  important 
role  in  this  experiment  and  therefore  it  has  not  been  considered  further 

As  far  as  boundary  reverberation  is  concerned,  it  was  possible  through  analysis  and 
interpretation  of  the  data  to  derive  the  conclusion  that  surface  reverberation  effects 
were  lower  than  one  might  have  expected.  It  was  also  concluded  that  the  experi¬ 
mental  set-up  was  not  sensitive  to  surface  scattering  and  therefore  it  was  essentially 
not  useable  for  estimating  surface  scattering  strength  values.  On  the  other  hand 
the  measured  low-frequency  reverberation  data  were  clearly  dominated  by  bottom 
scattering,  which  in  turn  affected  the  returns  in  the  various  ways  described  in  this 
section  Prior  to  the  description  however,  a  discussion  is  given  of  the  reasons  that 
led  to  the  rejection  of  the  surface-scattering  data 

S.l.  SURFACE  AND  BOTTOM  BACKSCATTERING  EFFECTS 

The  approach  taken  in  the  identification  and  separation  of  the  surface  and  bottom 
contributions  had  two  elements:  (a)  comparison  with,  and  analysis  of,  computer 
model  predictions,  and  (b)  clues  provided  by  linear  scattering  theory. 

(a)  The  computer  model  analysis  was  mainly  made  with  the  GENERIC'  Sonar 
Model  [9]  since  this  was  regarded  as  the  most  complete  model  available  for 
reverberation  predictions.  More  recent  programs,  such  as  the  NRL  Rever¬ 
beration  Model  are  possibly  more  tailored  to  the  analysis  of  the  data  that 
were  obtained,  but  they  were  unavailable  at  the  time  of  the  analysis. 

(b)  With  regard  to  the  search  for  clues  provided  by  linear  scattering  theory,  it 
was  decided  that  as  criteria  for  separating  the  two  processes  time  and  fre¬ 
quency  discriminants  could  be  used.  Since  long-duration  waveforms  were 
used  (minimum  pulse  width  1  s),  separation  in  time  was  really  only  appli¬ 
cable  before  reception  of  the  first  bottom  bounce.  After  that,  surface  and 
bottom  returns  overlapped.  Frequency  thus  seemed  possibly  to  be  the  best 
discriminant.  Since  the  conditions  of  the  experiment  were  such  that  the 
Rayleigh  parameter  was  less  than  unity,  the  theory  in  Subsect  3.4  indicated 
that  surface  scatter  ought  to  be  doppler  shifted  by  a  fixed  amount  dependent 
on  the  carrier,  whereas  bottom  scatter  would  have  to  be  at  zero  doppler  shift 
(assuming  towship  motion  compensation  to  be  applied) 

The  version  of  the  GENERIC  Sonar  Model  that  was  available  (version  C)  repre¬ 
sents  a  considerable  improvement  over  the  older  N1SSM  model,  but  nevertheless 
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retains  some  limitations  which  must  still  be  taken  into  account  when  comparing 
predictions  with  real  data.  The  most  relevant  are  the  following:  the  model  assu¬ 
mes  an  environment  independent  of  range  and  time,  that  is,  a  flat  bottom  and  a 
horizontally  stratified  ocean  (constant  sound  speed  profile  vs  range).  Sloping  of 
the  bottom,  which  has  been  recognized  as  an  important  factor  for  bottom  rever¬ 
beration  [11],  is  not  accounted  for.  Additionally,  GENERIC  does  not  adequately 
model  the  three-dimensional  beam  patterns  of  a  linear  array,  and  does  not  have 
bistatic4-e.  non-colocated,  source  and  receiver  modelling  capability.  Nevertheless, 
since  comparisons  have  been  purposely  made  with  appropriate  environments  such 
as  the  Abyssal  Plains,  the  model  indications  are  considered  acceptable.  An  example 
of  the  GENERIC  Sonar  Model  predictions  tailored  to  the  experimental  conditions, 
is  provided  at  Pig.  9,  in  which  the  reverberation  predictions  are  presented  for  the 
set-up  used  during  the  September  198b  Sea  Trial,  with  simulated  transmission  of 
an  8  s,  370  Hz  CW  pulse  in  sea-state  4  and  reception  at  the  broadside  beam. 

Inspection  of  the  results  obtained  with  signals  that  had  good  range  resolution  (tens 
of  meters)  shows  that  before  reception  of  the  first  bottom  bounce  the  data  in  all 
beams  appear  to  include  ambient  noise  alone,  i.e.  the  presence  of  surface  reverbe¬ 
ration  is  not  evident.  See  for  instance  data  at  the  initial  ranges  in  Fig.  13. 

After  the  onset  of  the  first  bottom  bounce,  reverberation  returns,  although  decaying 
as  a  function  of  distance,  were  up  to  50  dB  above  the  ambient  noise  level  at  the  initial 
ranges.  The  closer  that  beams  were  to  broadside  the  higher  the  received  levels.  The 
fact  that  these  beams  were  receiving  backscattered  energy  from  the  seafloor  more 
directly  (because  of  the  conical  symmetry  of  the  patterns)  is  an  indication  that 
bottom  reverberation  dominated  other  forms  of  scattering.  If  surface  scatter  were 
to  dominate,  the  intensities  would,  contrary  to  the  experimental  observations,  be 
increasing  with  increasing  angle  from  broadside,  particularly  in  the  convergence 
zones,  which  is  also  shown  by  model  predictions. 

Conclusions  drawn  from  the  received  power  levels  are,  however,  not  too  reliable — 
mainly  because  the  fathometer  effect  (to  be  discussed  later)  may  well  overwhelm  all 
other  processes,  but  also  because  the  bottom  and  surface  reverberation  levels  are 
not  too  different,  as  indicated  by  computer  predictions.  The  predictions  at  Fig.  9 
were  obtained  by  entering  scattering-strength  values  from  the  standard  formulae 
of  Chapman- Harris  and  McKenzie.  It  can  be  seen  that  although  bottom  returns 
almost  always  exceed  surface  returns,  except  in  the  convergence  tone,  the  difference 
is  consistently  marginal  (only  a  few  dB).  This  is  because  the  Chapman-Harris  for¬ 
mula,  unlike  linear  theory,  predicts  a  strong  dependence  of  5,  with  sea-state  (or 
wind  speed).  From  sea-state  2  to  sea-state  4  the  predicted  increase  of  the  scattering 
strength  at  medium  and  low  grazing  angles  wasis  approximately  25  dB. 

The  use  of  frequency  as  a  discriminant  relies  on  the  analysis  of  waveforms  with  high 
doppler  resolution.  The  spectral  analysis  of  long  CW  signals  does  not  indicate  the 
presence  of  doppler  shifted  peaks  in  the  received  spect--’  time  history  of  the  va- 
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rious  beams  which  would  he  predicted  if  surface  scattering  were  present.  Refer  for 
example  to  Figs.  22  to  25  These  are  contour  plots  of  ensemble-averaged  spectral 
time  histories  of  the  returns  following  an  R  s,  ;170  Hz  CW  transmission  The  figures 
show  data  for  two  beams,  broadside  and  rear  endfire,  and  for  two  different  sea- 
states:  sea-state  2  (average  wind  speed  less  than  7  kn)  and  sea-state  4  (wind  speed 
>  20  kn).  Although  the  locations  on  the  two  days  were  different,  it  was  essentially 
the  same  environment  (as  regards  seafloor  depth  and  structure)  and  velocity  pro¬ 
file  (summer;  Mediterranean;  deep-water).  It  is  apparent  from  the  broadside  data 
that  the  peaks  of  the  spectra  all  lined  up  at  zero  doppler  shift  and  that  there  is 
no  significant  secondary  lobes  at  the  frequency  shift  of  0.8  to  0.9  Hz,  as  postulated 
by  the  diffraction  grating  equation  (see  Eq.  (18))  At  endfire,  the  received  levels, 
initially  much  lower  than  the  corresponding  broadside  data,  show  shifted  spectral 
maxima  at  the  initial  ranges  which,  however,  soon  move  to  the  zero  frequency  shift 
position.  As  at  broadside,  there  is  also  at  endfire  a  lack  of  significant  secondary 
doppler  shifted  lobes.  In  all  the  figures  the  insets  reaffirm  these  findings  bv  showing 
the  spectrum  for  a  particular  range  at  which  a  possible  maximum  of  the  surface 
reverberation  could  have  been  observed,  as  predicted  by  computer  modelling.  This 
corresponds  to  the  first  convergenge  zone.  The  spectra  are  also,  in  this  instance, 
centered  at  zero  doppler  shift  and  are  also  quite  symmetrical,  with  no  pronounced 
side  peaks.  The  figures  thus  confirm  the  hypothesis  that  the  bulk  of  the  back- 
scattering  conies  from  the  seafloor.  Bottom  scattering  is  received  directly  at  the 
broadside  beams  and  is  received,  in  the  beams  towards  endfire,  initially  through  the 
sidelobes,  and  then  in  the  mainbeam  when  scattering  originates  from  that  portion 
of  the  seafloor  that  intersects  the  main  beam.  This  readily  explains  the  doppler 
shifted  initial  data,  which  are  due  to  the  lack  of  doppler  compensation  by  the  ODN 
algorithm  (see  the  discussion  at  Subsect.  2.4).  This  is  also  in  good  agreement  with 
model  predictions. 

Another  important  finding  that  points  to  a  surface  scatter  lower  than  predicted 
by  standard  formulae,  is  the  general  observation  that  the  returns  are  essentially 
independent  of  not  only  the  sea-state  (as  shown  in  the  previous  figures)  but  also 
of  the  sound  velocity  profile.  Figure  26  shows  returns  recorded  under  very  simi¬ 
lar  conditions  but  in  different  seasons  (winter  and  summer  profiles)  and  yet  they 
produced  very  similar  results.  Note  that  the  winter  conditions  ought  to  emphasize 
surface  reverberation  because  the  SSP  is  essentially  totally  upward  refracting. 

In  summary,  it  is  concluded  that  surface  reverberation  did  not  affect  the  pre¬ 
sented  results  and  that  scattering  from  the  seafloor  dominated  returns  However 
the  surface  still  had  a  profound  effect  on  the  received  process  by  modulating  the 
forward  and  backward  travelling  rays  associated  with  bottom  scattering  (Roderick 
and  Cron  [25]).  This  explains  the  high  tails  of  the  sidelobes  of  the  received  spec¬ 
tra.  These  are,  however,  centered  at  zero  doppler  shift  and  mostly  have  a  3-dB 
bandwidth  (after  matched  filtering)  within  the  expected  theoretical  limits 

Bottom  backscattering,  which  dominates  returns  of  low-frequency,  high-power,  long 
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duration  waveforms,  is  a  complicated,  non- planar  process.  It  appears,  however, 
to  be  quite  clearly  separated  into  two  distinct  processes,  in  agreement  with  the 
boundary  scattering  theory  previously  discussed.  These  are 

(a)  a  near- specular,  reflected  component  referred  to  as  the  fathometer  effect, 
and 

(b)  a  diffuse,  incoherent  component  representing  bottom  backscattering  coming 
from  progressively  smaller  grazing  angles. 

The  separation  of  the  two  contributing  processes  is  evident  only  if  the  transmitted 
waveform  has  enough  resolution  capability  in  time.  To  illustrate  this  for  a  10  Hz 
linear  PM  signal,  a  typical  example  of  the  two  separate  processes  is  presented  as 
Fig.  27a, b  (from  Ref.  [26]),  which  plots  the  ensemble  averaged  received  intensity 
time  history  recorded  at  the  broadside  beam  (90°)  and  at  beam  direction  128° 
(from  forward  endfire),  over  the  initial  20  km.  The  effects  associated  with  the 
contributing  processes  are  readily  evident  in  the  figure:  one  is  the  regularly  spaced 
peaks  whose  positions  are  beam  independent,  and  the  other  one  is  the  progressive 
decay  of  the  received  levels  as  a  function  of  time,  observable  in  between  the  peaks. 
For  the  128°  beam,  the  decay  of  the  received  power  levels  vs  time  associated  with 
the  seafloor  scattering  at  progressively  smaller  grazing  singles  only  begins  after  the 
backscattered  energy  is  received  into  the  main  beam  (the  peak  is  marked  by  an 
arrow  on  Fig.  27b). 


5.2.  THE  FATHOMETER  EFFECT 

The  fathometer  effect  arises  from  the  energy  that  is  radiated  by  the  source  in 
the  near-vertical  direction,  and  that  propagates  vertically  in  the  water  column, 
repetitively  bouncing  off  the  bottom  and  the  surface.  This  gives  rise  to  the  regularly 
spaced  peaks  (the  fathometer  peaks),  whose  range  separation  is  equal  to  the  water 
depth.  The  peaks  are  clearly  visible  in  Fig.  27a  and  in  the  FM  displays  of  all  the 
figures  referred  to  in  Sect.  4.  The  fathometer  peaks  are  obviously  more  evident  at 
the  broadside  beam,  but  are  sdso  present,  although  attenuated,  through  the  spatial 
sidelobes,  in  all  other  beams,  as  shown  at  Fig.  27b. 

The  fathometer  effect  essentially  involves  a  reflection  process  with  the  elementary 
areas  of  the  ensonified  surface  that  are  normal  to  the  incident  sound.  Then  to  a  first 
order  of  approximation  a  loss  factor  is  all  it  takes  to  describe  the  process.  Indeed, 
the  received  level  at  the  fathometer  peak  can  be  expressed  by 

P„  =  SL  -  TL„  -n£b  ,  (27) 

where  P„  is  the  received  power  level  at  the  nth  bottom  bounce  (dBrepPa);  SL  is 
the  source  level  (dBrepPa);  TL„  is  the  transmission  los-  (dB)  =  201og(2ndw)  for 
spherical  spreading:  d,  is  the  water  depth  (m);  Lb  is  the  bottom  loss  at  normal  in- 
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cidence  (dB).  It  has  been  assumed  that  the  surface  loss  and  volumetric  attenuations 
are  negligible,  since  frequencies  are  well  below  1  kHz, 

Thus  Lb  can  be  evaluated,  for  instance,  by 

A r  E(P„)  201og(~-j  +  U .  (2«) 

where  the  decay  between  successive  fathometer  peaks,  A p,  is  obtained  by  taking 
i  nst  mbli  averaged  consent! ive  fathometer  peaks 

Normal- incidence  bottom- reflect  ion  loss  estimates  made  in  this  way  indicate,  for 
the  measurement  area  in  the  Balearic  Basin,  values  of  5  to  8  dB  per  bounce  (at 
370  Hz),  which  are  in  good  agreement  with  acoustic  reflectivity  measurements  made 
in  the  same  areas  (Hastrup  and  Akal  [27]). 

5.3.  BOTTOM  BACKSCATTERING  US  GRAZING  ANGLE 

Since  bottom  reverberation  is  the  principal  contributor  to  the  received  reverbera¬ 
tion  levels,  use  can  be  made  of  the  received  intensity  time  history  to  determine 
bottom  scattering  strength  versus  grazing  angle.  For  that  purpose,  however,  it  is 
necessary  to  avoid  the  contaminations  caused  by  the  fathometer  effect.  Marandino 
and  Goldsberry  (26]  discuss  a  method  of  determining  bottom  scattering  strength 
versus  grazing  angle  by  taking  advantage  of  the  multiple  beamforming  capabilities 
of  the  towed  array  and  processor.  This  technique  offers  the  advantage  of  evaluating 
scattering  strength  at  local  maxima  of  the  ratio  of  reverberation  to  ambient  noise 
while  at  the  same  time  avoiding  fathometer  contamination.  The  method  requires 
the  use  of  waveforms  with  relatively  high  time  resolution,  capable  of  clearly  sepa¬ 
rating  surface  bounces  from  the  direct  returns  and  using  information  in  the  beams 
that  are  off  broadside.  The  reference  describes  this  technique  in  detail  and  gives  an 
example  of  the  evaluation  of  the  bottom  scattering  strength  versus  grazing  angle 
at  370  Hz  for  an  area  of  the  Balearic  Basin.  The  same  result,  expanded  to  include 
measurements  made  in  the  same  location  in  the  740  Hz  band,  is  given  at  Fig.  28. 
These  results  are  presented  here  in  support  of  conclusions  that  are  not  derived  from 
this  data  set  alone.  First,  the  general  behaviour  of  5b  versus  grazing  angle  indi¬ 
cates  again  that  for  seafloor  scattering  a  composite  roughness  model  is  applicable, 
with  a  strong  reflective  component  in  the  near-specular  region  and  a  diffuse  process 
for  medium  to  low  grazing  angles,  where  the  scattering  strength  does  not  seem  to 
depend  much  on  the  grazing  angle.  A  Lambert’s  rule  does  not  adequately  describe 
the  medium  and  low  grazing  angle  region  with  a  single  parameter.  A  best  estimate 
for  p  (see  Subsect.  3  5),  acceptable  down  to  about  20“,  is  -34  dB.  This  value, 
although  within  range  of  reported  values,  is  lower  than  expected  for  the  scattering 
from  the  Abyssal  Plains.  Finally,  the  data  indicate  that  the  scattering  strength 
in  the  sub-kilohertz  region  shows  no  frequency  dependence,  which  is  in  agreement 
with  the  results  of  Berkson  et  al.  [28]. 
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5.4.  COMPARISON  WITH  MODEL  PREDICTIONS 

As  previously  pointed  out,  the  GENERIC  Sonar  Model  [9j  was  considered  the 
most  appropriate  model  for  reverberation  predictions  in  relation  to  the  present 
experimental  set-up,  and  has  therefore  been  extensively  used  for  simulations. 

Of  more  value  than  the  exact  numerical  forecasting  of  reverberation,  which  is  bound 
to  be  highly  dependent  on  the  detailed  environmental  and  systematic  related  con¬ 
ditions,  is  the  model's  provision  of  a  key  to  the  parameters  which  are  relevant  in 
the  predictions. 

The  version  of  GENERIC  is  essentially  bi- dimensional  and  thus  does  not  allow 
precise  modelling  of  the  effect  of  a  three-dimensional  beam  pattern  on  the  returns 
Therefore  no  agreement  is  found  for  comparison  with  data  from  endfire  beams, 
although  surface  reverberation  is  possibly  correctly  modelled.  For  broadside  beam 
data  however,  where  a  two-dimensional  model  is  adequate,  the  agreement  is  gene¬ 
rally  more  satisfactory. 

Because  of  the  significance  of  the  scattering  from  the  seafloor,  modelling  is  use¬ 
ful  in  revealing  how  critical  the  correct  estimation  of  both  forward  (bottom  loss) 
and  backward  (scattering  strength)  parameters  is  for  obtaining  a  good  agreement 
between  measured  and  predicted  values.  Whenever  the  forward  and  backward  para¬ 
meters  have  been  completely  determined  the  agreement  is  particularly  satisfactory, 
as  shown  in  the  Figs.  29  and  30.  These  figures  show  GENERIC  predictions  su¬ 
perposed  on  actual  data,  for  reverberation  returns  measured  in  the  Balearic  Basin 
under  summer  and  winter  conditions,  and  for  two  waveforms  (8  s  CW  and  10  Hr 
LFM)  received  at  the  broadside  beam;  use  was  made  of  the  bottom  backscattering 
values  shown  at  Fig.  28. 

5.5.  LONG-RANGE  BACKSCATTERING 

Long-range  returns  are  dominated  by  directional  returns  from  highlights  of  the 
seafloor,  such  as  seamounts,  sea  pinnacles  and  islands,  coastlines  and  continental 
shelves.  This  dominance  is  seen  in  the  real-time  colour  displays  of  Sect.  4,  where 
the  Sardinian  and  Corsican  shelves  stand  out  clearly  for  measurements  made  in  the 
Balearic  Basin,  and  for  measurements  in  the  Tyrrhenian  Sea,  where  the  island  of 
Ustica,  the  Pliny  and  Vavilov  seamounts,  and  the  Sicilian  and  Neapolitan  coasts 
provide  large  returns.  The  directionality  is  more  visible  with  LFM  waveforms, 
which  have  high  range  resolution  capability,  but  is  equally  valid  for  CW  signals. 
Estimates  of  the  equivalent  target  strength  of  the  major  bathymetrical  features  can 
be  as  high  as  60  dB. 

Good  correlation  for  the  major  features  is  found  between  the  LFM  and  CW  re¬ 
turns,  with  the  latter  all  centered  at  zero  doppler  shift,  indicating  that  they  are 
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Attributable*  to  stationary  targets  All  doppler-shifted  returns  discernible  in  the 
CW  displays  can  easily  be  associated  with  major  features  in  other  directions  Since 
these  returns  are  picked  up  through  the  sidelobes,  the  beamforined  output  is  not 
only  attenuated  but  also  frequency  shifted,  as  discussed  in  Suhsect.  2.4  Large 
features,  with  a  high  equivalent  target  strength,  generate  echoes  50  dB  or  more 
above  the  ambient  noise  level,  and  since  practical  sidelobe  rejection  is  no  better 
than  30  dB,  on  beainfornung  these  returns  will  afreet  any  other  beams  with  levels 
of  15  dB  or  inure  above  ambient  noise.  All  returns  from  major  bathymetric  features 
have  a  complex  sidelobe  structure,  both  in  range  (or  time)  and  in  frequency;  this  is 
most  likely  a  consequence  of  the  multipath  associated  with  the  long-range  returns. 

Figures  31  to  34  show  various  examples  of  the  spectral  time  histories  of  returns 
from  transmissions  of  an  K  s,  370  Hz  signal,  along  beam  directions  with  various 
major  features  such  as  seamounts  and  pinnacles  (for  examples  the  Pliny  and  Vavilov 
seamounts  at  Figs.  31  and  32  respectively ),  islands  (Usticaat  Fig.  33)  and  coastlines 
or  continental  shelves  (western  coast  of  Corsica  at  Fig.  34).  These  figures  are  plots 
of  the  bathymetry  along  the  beam  pointing  direction,  all  with  a  common  range  scale 
adjacent  to  the  received  spectral  time  history  in  order  to  illustrate  the  correlation 
between  received  intensities  and  bottom  geometrical  features  (depth  and  sloping). 
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6.  Statistical  analysis  of  the  received  intensities 


«.l.  INTRODUCTION 

The  interest  in  a  higher-order  statistical  analysis  of  the  received  intensities  for 
an  activated  towed  array  is  twofold:  on  the  one  hand  it  helps  to  gain  a  better 
insight  into  long-range  propagation  and  scattering  effects  of  importance  in  the  sub¬ 
kilohertz  region;  on  the  other  hand  it  provides  a  useful  input  to  the  evaluation  of 
the  performance  of  detection  and  parameter  estimation  algorithms  suitable  for  the 
possible  operational  applications  of  the  activated  towed  array  concept. 

Reverberation  is  a  highly  non-stationary  process,  and  therefore  the  only  strictly  va¬ 
lid  statistics  are  those  that  apply  to  ensemble  averaging.  However  it  is  not  always 
possible,  for  practical  reasons,  to  collect  large  sets  of  independent  sample  transmis¬ 
sions  (pings)  for  ensemble  averaging.  Indeed,  in  the  experimental  set-up  the  number 
of  pings  over  wltich  one  estimates  ensemble  statistics  is  always  small  (a  maximum 
of  30),  rind  cannot  be  made  much  larger  if  homogeneity  among  data  sets  is  to  be 
preserved;  the  variation  among  data  sets  is  due  to  the  changing  of  the  environment 
caused  by  the  motion  of  the  acquisition  platform.  Additionally,  from  the  long-range 
sonar  detection  standpoint,  decisions  as  to  the  presence  or  absence  of  targets  must 
usually  be  made  on  the  basis  of  small-sample  sets.  Therefore  one  first  establishes 
the  relationship  between  the  small-sample  sets  and  their  ensemble  statistics  in  order 
to  derive  the  general  properties  of  the  medium  and  evaluate  detection  algorithms. 

The  analysis  carried  out  so  far  has  investigated  the  amplitude  statistics  of  the  re¬ 
ceived  intensities  after  beamforming,  matched  filtering  and  envelope  detection.  It 
is  therefore  essentially  a  narrowband  variance  analysis.  An  example  of  the  input  is 
given  at  Figs.  35  and  36.  The  waveforms  that  have  been  analysed  had  relatively 
high  range  resolution  (LFMs  or  short  CWs)  and  yielded  enough  independent  time 
samples  for  reliable  higher-order  statistics.  Figure  35  shows  a  typical  time  history 
of  the  received  intensities  for  one  single  ping  and  a  LFM  signal  of  8  s  duration  and 
10  Hz  bandwidth  (broadside  and  endfire  beam  data).  Figure  36  presents  the  same 
data  as  Fig.  35,  but  all  pings  for  the  same  run  have  been  ensemble  averaged.  The 
smoothing  effect  is  evident  although  the  number  of  pings  is  limited  (only  15).  En¬ 
semble  linear  averaging  is  applied  by  the  combining  and  averaging  of  the  magnitudes 
from  different  pings  sampled  at  the  same  time  instant  relative  to  their  individual 
transmission  instants.  The  standard  deviation  of  the  intensities  (expressed  in  dB) 
is  related  to  the  coefficient  of  variation  of  the  linear  averaging  and  is  indicated  in 
the  lower  half  of  each  plot  at  Fig.  36,  with  the  scale  on  the  right-hand  side.  Since 
the  tow  vessel  is  moving  during  the  run,  features  that  are  independent  of  actual 
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location  will  be  properly  averaged  whereas  those  that  are  not  independent,  such  as 
returns  from  objects  whose  relative  geometry  changes,  will  he  smeared. 

In  Figs.  35  and  36  it  is  possible  to  distinguish  three  discrete  regions  dominated 
by  three  distinct  processes.  The  first  region,  which  extends  up  to  about  40  kin, 
is  dominated  by  reverberation,  in  particular  bottom  reverberation  as  discussed  in 
the  previous  sections.  The  second  region,  which  extends  beyond  the  first  one  up  to 
about  90  km,  is  dominated  by  large  returns  from  an  extended  target,  here  part  of 
the  Continental  shelf  off  the  western  coast  of  Corsica.  The  third,  beyond  90  km, 
appears  to  be  dominated  by  ambient  noise  alone.  This  clear-cut  separation  made 
this  kind  of  data  set  convenient  for  the  present  statistical  analysis  and  has  been 
used  in  the  following  sections  to  illustrate  the  method  of  application. 

The  methodology  is  based  on  the  extensive  use  of  non-parametric  hypothesis  te¬ 
sting,  particularly  in  regard  to  stationarity  and  the  distribution  types  of  the  received 
intensity  time  series,  since  these  are  the  most  relevant  statistics  for  sonar  detection. 
The  advantage  of  using  non-parametric  statistics  is  that  they  are  ‘distribution- 
free’,  i.e.  they  do  not  need  an  underlying  assumption  (usually  of  normality)  on 
the  population.  This  is  important  when  either  the  number  of  samples  is  small  or 
when  the  data  are  by  their  nature  non-normal  or  non-stationary,  as  is  the  case  in 
reverberation-dominated  processes.  Note  too  that  in  order  to  apply  the  statistical 
tests  it  is  necessary  to  operate  on  independent  samples.  Since  the  processed  data 
base  is  oversampled,  a  decimation-in-tiine  is  required.  Decimation  by  a  factor  of  two 
usually  gives  enough  de-correlation  to  provide  the  required  independence  between 
successive  samples. 

In  the  following  section,  non-parametric  tests  are  introduced  for  the  analysis  of  data 
limited  by  (stationary)  ambient  noise.  Since  reverberation-dominated  stochastic 
processes  can  usually  (with  care)  be  reduced  to  stationary  processes,  as  it  will 
be  shown  later,  the  same  statistical  toois  can  be  applied  to  the  analysis  of  the 
fluctuations  process  obtained  by  the  decomposition  of  the  original  process. 

6.2.  ANALYSIS  OF  AMBIENT  NOISE  LIMITED  DATA 

Although  ambient  noise  intensities  below  1  kHz  are  subject  to  long  term  variations 
and  are  significantly  affected  by  the  level  of  shipping  noise,  it  is  accepted  that  over 
limited  time  intervals  (a  few  minutes)  stationarity  can  usually  be  assumed,  see  for 
instance  [29].  Examination  of  the  data  set  under  investigation  (shown  at  Figs.  35 
and  36)  revealed  it  to  be  well  behaved  and  accordingly  provided  a  good  basis  as  a 
reference  for  the  application  of  the  techniques.  The  good  behaviour  was  attributed 
to  there  being  no  nearby  or  remote  shipping  affecting  returns  over  all  pings,  and  no 
contamination  of  beams  by  towship  noise  other  than  in  the  forward  endfire  beam. 

The  grand  average  of  the  intensities  received  over  those  ranges  dominated  by  ambient 
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noise  is  consistent  with  that  measured  independently  on  a  broader  band  with  a 
single  hydrophone,  and  is  in  good  agreement  with  the  Wenz/Knudsen  curves  for 
the  wind  and  sea-state  conditions  at  the  time  of  the  run.  The  average  value  is 
46.8  dBre/iPa  (see  Table  6),  and  by  adding  25  dB  for  array  and  processing  gain  a 
noise  level  of  71.8  dB  re/iPa/Hz  is  obtained. 

But  the  point  of  interest  is  the  fact  that  in  general  the  ensemble  average  of  the  stan¬ 
dard  deviation  of  the  dB  averages  over  time,  is  5.56  dB.  This  result  is  in  agrement 
with  standard  theories  of  long-range  ocean  acoustic  propagation  in  a  SOFAR-like 
channel  (see  Frisk  [301).  Notably,  Dyer's  theory  [31j  shows  that  if  scattering  can  be 
considered  negligible  with  respect  to  propagation  effects,  multipath  can  be  model¬ 
led  as  an  incoherent  superposition  of  many  (although  as  few  as  4  may  work  as  well) 
random  vectors  of  approximately  identical  energy  associated  with  each  path.  On 
converting  the  received  intensities  to  a  logarithmic  scale,  therefore,  their  dB  mean 
has  a  standard  deviation  with  an  expected  value  of  5.57  dB. 

6.3.  STATION ARITY  TESTS 

The  indication  that  it  is  possible  to  interchange  time  and  ping  (or  ensemble)  avera¬ 
ges  for  first-order  statistics  points  to  a  form  of  a  weak  stationarity.  This  was  tested 
by  means  of  the  Mann- Whitney  V  Test,  which  is  generally  regarded  as  one  of  the 
most  powerful  non-par ametric  tests  for  the  case  of  two  independent  samples  (Sie¬ 
gel  [32])  and  avoids  many  of  the  weaknesses  of  the  Parametric  t  Test.  It  is  used  for 
testing  whether  two  groups  have  been  drawn  from  the  same  population,  through 
checking  on  the  population  medians,  which  for  skewed  distributions  is  usually  a 
more  reliable  indication  than  checking  the  mean. 

In  the  application  described  here  the  Mann-Whitney  U  Test  was  used  to  determine 
whether  the  received  intensities  over  the  assumed  noise-limited  region,  formed  a 
locally  stationary  sequence  on  a  single  ping  basis.  Since  in  sonar  detection  the 
interest  is  in  maintaining  the  properties  over  as  small  an  interval  as  is  feasible,  the 
methodoly  that  was  applied  was  to  divide  the  total  data  set  into  smaller  contiguous 
subsets  and  run  the  Mann-Whitney  U  Test  on  each  group  by  taking  each  subset 
pair  once  in  combination  with  all  the  others.  As  is  well  known,  the  procedure 
for  a  non-parametric  test  is  based  on  standard  hypothesis  testing  and  requires  the 
setting  of  a  threshold  g  (the  significance  level)  against  which  results  are  compared 
in  order  to  decide  whether  the  hypothesis  must  be  rejected  or  not.  In  the  present 
case  the  hypothesis  to  be  tested  was  the  stationarity  of  the  time  series  and  various 
significance  levels  were  tried  out.  However,  results  are  given  only  for  the  customary 
5%  level,  which  was  quite  representative 

Selecting  a  significance  level  g  means  that  if  the  test  can  be  run  for  many  times 
on  independent  samples  and  the  hypothesis  is  true,  there  should  be  approximately 
100p%  rejections.  In  the  stationarity  tests  the  total  number  of  independent  data 
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per  ping  was  fixed,  and  so  a  trade-off  resulted  between  the  subset  (or  record)  size 
and  the  number  of  independent  tests.  Tests  were  then  run  for  three  record  lengths 
on  each  ping.  The  results  are  shown  in  Table  4,  where  the  percentage  rejection 
obtained  by  running  the  tests  at  the  5%  significance  level  is  indicated.  The  results 
indicate  that  stationarity  cannot  be  rejected.  In  fact,  the  selected  testing  criterion 
is  rather  restrictive  in  that  it  is  equally  affected  by  records  which  are  close  t»r 
that  are  widely  separated  in  time.  The  test  is  quite  sensitive  even  to  relatively 
small  variations  of  the  rms  intensity  fluctuations,  and  therefore  in  the  presence  of 
reverberation  or  extended  target  returns  the  percentage  rejection  is  very  high. 


ii.4.  PROBABILITY  DISTRIBUTION  ESTIMATION 

Due  to  there  not  being  enough  samples  for  ensemble  statistics  probability  densities 
and  distributions  have  been  estimated  here  only  for  independent  turn  samples. 
However,  when  the  time  series  is  stationary  these  statistics  are  also  the  ensemble 
statistics.  The  evaluation  of  the  distributions  was  also  based  on  non-parametric 
statistical  methods  instead  of  the  plotting  of  histograms  or  the  tabling  of  percent  iles. 
The  method  consisted  of  determining  whether  at  a  preset  significance  level  a  given, 
theoretical  distribution  fitted  the  data,  or  out  of  a  number  of  distributions  one 
fitted  the  data  best.  The  Kolmogorov-Smirnov  (KS)  goodness  of  fit  test  [32]  which 
was  used  for  this  purpose  is  considered  well  suited. 

The  theoretical  densities  that  were  taken  into  account  for  the  received  intensities 
(magnitude  data)  were 

(i)  The  Rayleigh  or  2-degrees-of- freedom  \2  density  (this  is  the  negative  expo¬ 
nential  density  if  the  data  are  expressed  as  power  values,  and  is  the  density 
of  the  envelope  of  a  zero-mean  gaussian  input  sequence). 

(ii)  The  log-normal  density,  which  does  not  arise  from  physical  considerations, 
but  is  consistent  with  letting  the  parameters  of  the  sonar  equation  to  be 
gaussian  random  variables. 

(iii)  The  Weibull  density,  Ref.  [33],  which  is  a  two-parameter  distribution  that  is 
sometimes  used  to  describe  radar  land  clutter.  It  is  a  distribution  that  can 
be  made  to  fit  either  a  Rayleigh  or  a  log-normal  density  or,  heuristicaliy, 
distributions  in  between. 

Other  distributions  can  be  considered,  such  as  that  arising  front  the  gamma  den¬ 
sity  and,  notably,  the  Rice  distribution  which  possibly  is  the  most  complete  if 
associated  with  forms  of  coherent  scattering  or  dominating  path.  However  the  Rice 
distribution  is  pretty  unwieldy  and  (Urick  [34])  it  requires  the  estimate  of  an  extra 
parameter,  which  ruled  it  out  for  this  study. 

The  criterion  in  the  application  of  the  KS  test  is  the  same  as  the  general  one  used  in 
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TABLE  4 

Results  of  the  Maim-Whitney  U  test  for  stationarity 


Stationarity 

( Mann-Whi tney  U)  tests 

Day :  6  Sept  85 

-  Signal:  Ml  -  Broadside 

Beam 

Range  min  89.96  Range  max 

156.26 

110  data/record 

4  records 

6  MW  U  tests 

Ping  ft  % 

rejection  (5%  level) 

1 

0.00 

2 

0.00 

3 

0.00 

4 

33.33 

5 

16.67 

6 

0.00 

7 

0.00 

8 

0.00 

9 

0.00 

10 

0.00 

11 

16.67 

12 

0.00 

13 

0.00 

14 

50.00 

15 

0.00 

avrg  %  rejection  at  5»  level  -  7.78 

55  data/record 

8  records 

28  MW  U  tests 

Ping  ft  % 

rejection  (5%  level) 

1 

3.57 

2 

7.14 

3 

3.57 

4 

7.14 

5 

14.29 

6 

0.00 

7 

0.00 

8 

0.00 

9 

0.00 

10 

0.00 

11 

0.00 

12 

3.57 

13 

17.86 

14 

7.14 

15 

0.00 

avrg  %  rejection  at  5%  level  -  4.29 

36  data/record 

12  records 

66  MW  U  tests 

Ping  ft  % 

rejection  (5%  level) 

1 

7.58 

2 

0.00 

3 

1.52 

4 

4.55 

5 

6.06 

6 

1.52 

7 

4.55 

8 

0.00 

9 

4.55 

10 

0.00 

11 

7.58 

12 

6.06 

/ 

13 

10.61 

14 

12.12 

15 

6.06 

avrg  %  rejection  at  5%  level  -  4.85 
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the  non-parametric  tests,  namely  whether  a  set  of  independent  samples  conies  from 
a  specified  theoretical  distribution.  The  KS  test  determines  the  maximum  difference 
(D)  between  the  theoretical  and  the  experimental  cumulative  distributions.  The 
sampling  distribution  of  D  is  known  under  the  hypothesis  of  acceptance,  and  thus 
the  estimated  value  can  then  he  compared  with  a  threshold  to  determine  whether 
to  reject  the  hypothesis  at  a  specified  significance  level. 

The  test  can  be  applied  on  each  ping  for  a  single  set  of  all  the  available  independent 
data  points.  The  results  for  the  data  set  presented  in  Table  4  are  given  along 
with  first  and  second  order  statistics  at  Table  5.  The  rms  intensities  are  definitely 
Rayleigh  distributed,  with  an  average  standard  deviation  of  the  dB  mean  of  5.56  dB. 
This  result  has  a  general  vlaidity.  Tests  performed  on  different  beams  and  runs  show 
that  if  the  set.  passes  the  stationarity  test  then  the  data  obey  Rayleigh  statistics. 
Note  also  that  such  tests  are  only  really  applicable  to  single-ping  data:  the  ensemble 
averaged  results  at  Fig.  36  (only  data  points  beyond  90  km)  have  a  very  small 
variance  (i.e.  have  a  dominant  deterministic  component),  and  thus  the  test  loses 
significance. 

The  experimental  set  of  stationary,  noise-limited  data  is  plotted  in  Fig.  37a  as  a 
function  of  the  normalized  amplitude;  Fig.  37b  shows  the  corresponding  cumulative 
distribution  function  versus  the  three  theoretical  distributions  introduced  previou¬ 
sly.  The  horizontal  axis  in  the  figure  is  the  amplitude  normalized  to  the  experimental 
median,  and  expressed  in  dB  to  facilitate  comparison  with  the  other  curves;  and 
the  vertical  axis  is  the  quantity  1  -  P(x),  in  log  scale,  where  P(z)  is  the  cumulative 
distribution  function.  This  rendition  allows  a  readier  appreciation  of  the  behaviour 
of  the  experimental  data  at  the  tails  of  the  distribution,  where  deviations  from  theo¬ 
retical  have  a  stronger  effect  on  the  probability  of  false  alarm.  An  optimum  linear 
processor  is  designed  to  operate  with  Rayleigh-distributed  envelopes  and  therefore 
higher  tails  (such  as  those  associated  with  a  log-normal  distribution)  cause  higher 
false  alarm  rates,  and  thus  worse  performance.  In  the  present  instance  it  can  be 
seen  that  the  data  closely  obey  Rayleigh  statistics.  This  being  the  case,  the  Weibull 
distribution  essentially  overlaps  the  Rayleigh  distribution. 


6.5.  REVERBERATION  STATISTICS 

When  the  process  includes  some  form  of  scattering  or  has  varying  rms  levels,  its 
statistics  are  highly  non-stationary  and  the  distributions  depart  from  the  Rayleigh 
model  and  tend  towards  a  log-normal  or  Weibull  distribution. 

Figure  38b  for  instance,  repeats  the  cumulative  distribution  function  versus  norma¬ 
lized  logarithmic  amplitude  for  the  data  set  shown  in  Fig.  38a — which  is  dominated 
by  reverberation  returns  (intensities  received  over  the  first  40  Km  of  range).  The 
large  departure  from  the  theoretical  Rayleigh  model,  and  the  agreement  over  the 
tails  with  a  Weibull  model  close  to  the  log-normal  case,  are  readily  evident.  Inte- 
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TABLE  5 

Kolmogorov-Smitnov  test  and  statistics  on  noise  limited  data 


Day:  6  Sept  1985 

-  Signal:  Ml  -  Broadside  Beam 

Range  rain  -  89. 

96  km,  Range  max  - 

156.26  km 

K  -  s  test  of  prob  density  type  - 

Rayleigh 

on  442  samples  per  ping 

Ping  • 

1 

dB  Mean  &  std  - 

46.4 

5.69 

Linear  mean  &  std  - 

251 . 

143, 

passed  K-s  Test  at 

2%  significance 

7? 

YE 

Passed  K-S  Test  at 

5%  significance 

?? 

YE 

Passed  K-S  Test  at 

10%  significance 

?? 

NO 

Ping  t 

2 

dB  Mean  &  std 

46.9 

5.60 

Linear  mean  &  std  - 

264. 

138. 

Passed  K-S  Test  at 

2\  significance 

?? 

YE 

Passed  K-S  Test  at 

5%  significance 

?? 

YE 

Passed  K-S  Test  at 

10%  significance 

?? 

YE 

Ping  1 

3 

dB  Mean  a  std  -  46.2 

5.91 

Linear  mean  &  std  -  245. 

130. 

Passed  K-S  Test  at 

21  significance 

?? 

YE 

passed  K-S  Test  at 

5%  significance 

?? 

YE 

Passed  K-S  Test  at 

10%  significance 

?? 

YE 

Ping  1 

13 

dB  Mean  &  std  -  46.6 

5.60 

Linear  mean  &  std  *  255. 

139. 

Passed  K-S  Test  at 

2%  significance 

?? 

YE 

Passed  K-S  Test  at 

5%  significance 

?? 

YE 

Passed  K-S  Test  at 

10%  significance 

?? 

YE 

Ping  t 

14 

dB  Mean  &  std 

.  47.3 

5.32 

Linear  mean  &  std  -  271. 

134. 

Passed  K-S  Test  at 

2%  significance 

?? 

YE 

Passed  K-S  Test  at 

5%  significance 

?? 

YE 

passed  K-S  Test  at 

10%  significance 

?? 

YE 

Ping  1 

15 

dB  Mean  &  std 

.  46.6 

5.15 

Linear  mean  &  std 

-  248 . 

125. 

Passed  K-S  Test  at 

2%  significance 

?? 

YE 

Passed  K-S  Test  at 

5%  significance 

7  7 

YE 

passed  K-S  Test  at 

10%  significance 

?? 

YE 

K-S  TEST  GLOBAL 

STATISTICS  on  15  pings 

avrg  dB  mean  - 

46.8  its  std 

0.357 

avrg  lin  mean  - 

259.  its  std 

11.0 

avrg  dB  std  -  5. 

56 

%  rejection  at  2% 

level  -  0.00 

%  rejection  at  5% 

level  -  0.00 

%  rejection  at  10% 

level  -  6.67 
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restingly,  although  the  Weibull  distribution  tils  the  data  better  at  the  higher  tails, 
the  KS  test  actually  indicates  that  the  data  adhere  more  closely  to  log-normal  sta¬ 
tutes.  This  is  because  the  test  considers  all  data  points,  whereas  the  (distorted) 
image  expands  a  relatively  small  but  important  percentage  of  data 

hollowing  the  approach  of  Ol’shevskii  |f>),  and  Plenums,  Shooter  and  Middleton  j;i.r>i, 
the  statistical  properties  of  reverberation  can  be  analysed  even  in  the  absence  of 
experimental  ensembles,  by  investigating  the  reverberation  time  covariance  and,  in 
particular,  reducing  the  process  to  a  (weak)  stationary  one  wherever  possible. 


The  covariance  of  a  stochastic  process  «(f),  is  given  by 

f, (/,,/,(  -  E{s{ f,  !-s( f 2 ) >  ,  (20) 

where  E{.  .)  denotes  ensemble  averaging.  For  a  weakly  stationary  process  is  a 
function  only  of  r  —  t\  -  f2.  If  the  stochastic  process  s (f)  can  be  represented  by 
the  product  of  two  independent  functions 

«(f)  =  r(f)tf(f)  (30) 

with  r(f)  a  rapidly  fluctuating,  stationary  process,  with  covariance  C,(r),  and  q{t) 
a  deterministic  function  representing  the  trend  or  slow  rms  variation,  it  follows  that 

C.(ti,t2)  =  <?(ft)<f(f2)G(r).  (31) 

And  then  if  the  function  q(t\  varies  slowly  over  the  averaging  interval  (which  occurs 
when  q(ti )  ^  q(t})) —  the  covariance  s(f )  is  a  process  reduceable  to  stationarity,  and 
the  decomposition  at  Eq.  (30)  is  applicable,  and  the  statistical  properties  can  then 
be  derived  by  time  averages.  The  local  covariance  stationarity  therefore  implies  the 
following: 

(i)  The  geometries,  and  beam  patterns  of  the  experiment,  as  well  as  the  scatte¬ 
ring  mechanisms  vary  slowly  over  the  time  intervals  of  averaging. 

(ii)  The  averaging  time  windows  must  be  appropriately  chosen. 


If  reverberation  can  be  assumed  to  be  locally  stationary,  then  a  technique  known  as 
logarithmic  normalization  (van  Schooneveld  (36) )  can  be  applied  to  decompose  the 
time  series  into  a  rapidly  fluctuating  and  stationary  process  and  a  slowly  varying 
function  called  the  trend,  which  represents  the  deterministic  contributions  of  the 
scattering  experiment.  By  taking  the  log  of  S(t )  as  expressed  by  Eq  (30)  it  is 


5(f)  =  R{t)  +  Q{t), 

where 

5(f)  =  10  log s(f) , 
R(t)=  10  log  r(f), 
<?(f)  =  101og?(f). 


(33) 
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The  logarithmic  conversion  transforms  the  normalization  into  a  linear  operation 
and  can  then  be  handled  through  linear  filtering  methods,  such  as  Fast  Fourier 
Transform  algorithms. 

Theoretically  an  analytic  description  of  the  deterministic  component,  such  as  the 
correlation  structure  or  its  spectrum,  should  be  available.  By  subtracting  the  spec¬ 
trum  from  the  frequency  representation  of  the  original  process  it  should  then  be 
possible  to  derive  the  properties  of  the  stationary  component.  This  is  seldom  fea¬ 
sible  in  practice;  however,  it  is  important  to  note  that  for  Eq.  (30)  to  hold  there 
must  be  an  appreciable  difference  between  the  correlation  characteristics,  and  thus 
between  the  spectral  behaviour  of  the  two  processes  (conditions  (i)  and  ( ii ) ) .  The¬ 
refore  although  the  analytic  structure  of  Q(t  )  may  not  be  known,  it  must  however 
be  confined  to  the  low-frequency  part  of  the  spectrum;  in  other  words  it  must  be 
a  slowly  varying  process.  If  that  is  the  case,  a  sub-optimum  technique  which  is 
simple  yet  quite  robust  can  be  used:  it  involves  choosing  a  threshold  wo  by  an 
iterative  process  which  leads  to  the  high-frequency  part  of  the  spectrum  passing 
stationarity  tests.  Having  set  such  a  threshold,  it  is  then  possible  to  apply  ideal 
filtering  by  removing  (through  setting  to  zero)  either  the  low  or  the  high  portion  of 
the  spectrum.  The  result  is  a  simplified  representation  of  the  type 

^5(w)  =  r/f(u>)  +  /-(?(w),  (33) 

where  the  prefix  T  represents  the  power  spectrum,  with 

J-  R(w  <  u>o)  =  0  , 


and 


3~Q(u>  >  u>o )  =  0 


(34) 


In  other  words,  instead  of  looking  for  the  optimum  transfer  function  to  decompose 
the  process,  a  simpler  spectrum  truncation  is  applied.  Note  that  this  is  a  batch 
processing  technique,  since  it  requires  having  the  whole  time  series  available. 

The  power  spectra  of  two  typical  sequences  of  received  log  intensities,  one  domi¬ 
nated  by  reverberation  and  the  other  dominated  by  ambient  noise,  are  shown  in 
Fig.  39.  Figure  39a  is  related  to  the  data  at  Fig.  37a,  ard  Fig.  39b  is  related  to 
the  data  at  Fig.  38.  The  frequency  range  of  the  transform  is  based  on  the  data 
sample  rate,  here  5  Hz  for  0.1  s  data  rate.  The  Fourier  transform  is  applied  by 
means  of  the  FFT  algorithm,  and  therefore  the  discontinuities  at  the  edges  of  the 
sample  introduce  some  noise  into  the  high-frequency  pert  of  the  spectrum.  When 
the  process  is  stationary  the  spectrum  shows  no  relevant  peaks  or  slopes.  When, 
on  the  contrary,  data  include  trend  variations  or  regular  patterns  coming  from 
gross  features,  the  power  spectrum  shows  the  presence  of  slopes  and  peaks  that 


Saclantcen  SR- 112 


affect  the  low-frequency  region.  It  is  in  general  visually  easy  to  distinguish  between 
low-  and  high-frequency  components  coming  from  slow  and  rapidly  varying  effects 
respectively. 

Bv  taking  the  inverse  Fourier  transform  it  is  thus  possible  to  extract  both  the  trend 
and  a  process  that  represents  the  rapid  fluctuation  around  the  trend  In  tile  log 
scale  this  last  process  is,  by  definition,  a  zero  mean  process;  its  statistics  are  of 
interest  to  the  sonar  detection  case.  The  linear  process,  which  has  a  normalized 
mean,  was  again  tested  with  the  KS  test  of  goodness  of  fit  to  determine  the  appli¬ 
cable  distribution  type.  The  general  result  is  that  even  in  the  presence  of  strong 
reverberation  a  well  chosen  cutoff  frequency  yields  data  with  a  strong  tendency 
to  be  Rayleigh  distributed.  Note  that  since  the  Rayleigh  distribution  is  a  single 
parameter  distribution,  and  since  for  the  normalized  process  the  median  is  unity 
by  definition,  the  variance  is  constrained.  An  example  of  the  fluctuations  process 
obtained  by  removing  the  trend  from  the  data  at  Fig.  38  is  shown  in  Fig.  40,  which 
also  presents  the  cumulative  distribution,  and  this  is  seen  to  closely  fit  the  Rayleigh 
model. 
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7.  Conclusions 

This  report  summarizes  results  of  reverberation  measurements  made  with  a  low- 
frequency  activated  towed  array  and  the  findings  of  a  preliminary  analysis  of  the 
data  collected. 

The  experimental  data  pertaining  to  the  reverberation  observed  when  operating  in 
the  Mediterranean  Sea  and  in  the  middle  of  the  sound  channel,  at  moderately  high 
power,  has  been  interpreted  and  provides  the  following  conclusions: 

(a)  Volume  reverberation  is  not  significant  for  the  frequencies  and  experimental 
set-up  utilized  in  this  locality. 

(b)  Surface  scattering  was  lower  than  predicted  by  the  standard  models,  and 
was  always  overshadowed  by  bottom  scattering. 

(c)  Bottom  reflections  in  the  vertical  water  column  (the  fathometer  returns)  and 
bottom  backscattering  dominated  the  received  intensity  time  series. 

(d)  Agreement  with  model  predictions  was  generally  good  whenever  the  model 
input  parameters  for  the  systemat-related  and  environmental  conditions  were 
accurately  known. 

(e)  Bottom  scattering  strength  was  reasonably  independent  of  frequency.  Spec¬ 
tra  with  little  or  no  doppler  shift  and  spread  were  observed,  but  with  high 
tails  in  the  sidelobes.  Long-range  returns  (beyond  about  20  n.mi)  was  usually 
dominated  by  directional  returns  from  extended  targets  such  as  seamounts 
and  pinnacles,  bottom  slopes  and  continental  shelves,  islands  and  coastlines. 

A  statistical  analysis  was  made  of  the  received  intensities  generated  by  transmitting 
frequency  modulated  and  relatively  short  (Is)  continuous  wave  signals.  Higher- 
order  statistics  were  investigated  by  extensive  use  of  non-parametric  hypothesis 
testing.  In  general,  it  turns  out  that  returns  dominated  by  ambient  noise,  if  not 
contaminated  by  towship  or  nearby  shipping  noise,  were  usually  stationary  and 
adhered  to  Rayleigh  statistics  to  various  levels  of  confidence.  Received  time  series 
dominated  by  reverberation  are  always  highly  non  stationary,  with  log-normal  or 
Weibull  amplitude  statistics  m  time,  however  reverberation  can  be  assumed  to  be 
locally  covariance  stationary.  By  applying  logarithmic  normalization  techniques  it 
is  thus  possible  to  separate  a  slowly  varying,  non-stationary  mean  from  a  rapidly- 
fluctuating  and  mostly  stationary  series  which  in  general  obeys  Rayleigh  statistics. 
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Fig  1.  Experimental  set-up  for  reverberation  measurements  with  an  activated 
towed  array. 


Fig.  2.  Geographical  areas  of  SACLANTCEN  sea  trials  for  deep-water  rever¬ 
beration  measurements. 
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Fig  I  Geophysical  and  acoustic  parameters  of  a  sample  of  seabed  in  the 
Balearic  Basin  (from  ■J'l 
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Fig.  4.  Typical  sound  speed  profiles  (early  April  and  early  September  data). 
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Fig  5  Shipborne  data -collection  system;  simplified  block  diagram 
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Fig.  6.  Filter  responses  for  370  liv.  m  (/uisiuon 
(a)  first  st Age  interpolator;  (b)  first  anti-aliasing  (in¬ 
terpolator);  (c)  second  anti  alining  f.-W'J. 
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Fig.  9.  GENERIC  Sonar  Mode I  predictions  (with  defauJt  values)  of  reverbe¬ 
ration  at  low-frequency,  simulating  this  experimental  set-up. 
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Fig.  10.  Surface  scattering  geometry  and  the  diffraction  grating  formula. 
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WAVE  NUMBER  (SURFACE  WAVES) 

Fig  11  Pierson-Mosiowily  sea  surface  wave  frequency  spectrum- 


Grazing  Angle  (deg)  Grazing  Angle  (deg) 

Fig.  12.  Composite-roughness  scattering,  (a)  surface  baciscattermg  strength 
(Chapman- Harris  formula);  (b)  bottom  baciscattering  strength  (McKenzie 
formula) 
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Fig.  13  Real-time  multi-ping  colour  display  of  received,  matched -filtered  intensities  as  a 
function  of  be  am  and  range  (refer  to  text  and  Table  3):  FM  display  only,  signal  C’F,  day  10, 
Apr.  '84  Data  in  dflrepPa. 
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Fig  IS.  Real-time  multi-ping  colour  display  of  received,  matched- filtered  intensities  as  a 
function  of  beam,  range  and  frequency  (refer  to  text  and  Table  3):  display  of  FM  and  CIV 
data ,  signal  Ml,  day  24,  Sept  '84.  Data  in  dBrtpPa 
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Fig.  17.  Real-time  mulli-ping  colour  display  of  received,  mat ched-filtered  intensities  as  a 
function  of  beam,  range  and  frequency  (refer  to  text  and  Table  3):  display  of  FM  and  C’VV 
data,  signal  M7,  day  4,  Sept  '85.  Data  in  dBrtpPa. 
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Fig  19  Real-time  multi-ping  colour  display  of  received,  matched- filtered  intensities  as  a 
function  of  beam,  range  and  frequency  (refer  to  text  and  Table  3):  display  of  FM  and  (’IV 
data,  signal  M7;  day  6,  Sept  '85-  Data  in  dBftpPa. 
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Fig  22.  Spectral  time  history  of  received  intensities :  8  s,  370  Hz  CW  tran¬ 
smission ;  broadside  beam .  sea-state  2. 
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Fig.  23.  Spectral  time  history  of  received  intensities:  8  s,  370  Hz  CW  tran¬ 
smission,  end  fire  be  am;  sea-state  2. 
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Fig.  24.  Spectral  time  history  of  received  intensities:  8  s,  370  Hz  OVV  tran¬ 
smission,  broadside  beam;  sea-state  4. 
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)an*>on  of  received  intensities  under  winter  (curves  IV, 
■>  S)  conditions:  6  s,  370  Hz  CIV'  pulse,  zero-th  frequenc) 
broadside  data:  (b)  end  fire  data. 
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Fig.  26.  Comp 
Summer  (curvet 
first  60  km  (a) 


SaC'LANTCEN  SR-U2 


SUMMER  85  BALEARIC  BASIN,  370HZ  CW  8SEC  .  17KN  HIND. 
2 1 3DB  DIPOLE  IX.  BROADSIDE  RX .  MY  DRTAi  80TIOS.BBS. 


a) 


SUMMER  OS  BALEARIC  BASIN.  FM  I DHZ :8SI SIM  I  -  I7KN  WIND 
370HZ  213D3  DIPOLE  IX.  BROADSIDE  RX.  MY  8OTL0S,  BBS  MY  DATA. 

b) 


TIME  IN  S 

Fig.  29.  Comparison  of  measured  reverberation  returns  -s  predictions  (with 
GENERIC  Model)  over  the  first  SO  im  of  range;  broadside  data;  summer 
environment  (a)  8  s,  CW  pu/se  (CW  data  from  Ml :249);  (b)  10-Ha  LFM 
pu/s<  (FM  data  from  M  1:249). 
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WINTER  04  BALEARIC  BR5IN.  370HZ  CW  85EC ,  22KN  HIND. 
209DB  OMNI  TX .  BROADSIDE  RX.  MY  DATA*  BOTL0S  .BBS 


a) 


WINTER  84  BALEARIC  BASIN.  FM  l OHZ  :8SEC f 51 M  1  .  20KN  WIND 
370HZ  209  DB  OMNI  TX.  BROADSIDE  RX.  MY  DATA  i  BOTL0S.  BBS. 


Fig.  30.  Comparison  of  measured  reverberation  returns  vs  predictions  (with 
GENERIC  Model)  over  the  first  60  km  of  range:  broadside  data;  winter  en¬ 
vironment  (a)  8  s,  CW  pulse  (CW  data  from  Ml:  101);  (b)  10-Hz  LF M  pulse 
(FM  data  from  AfJ.JOJ). 
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Fig.  31.  Spectral  time  history  and 
bathymetry  in  selected  beam  pointing 
directions  (bottom  highlight  here  is 
Pliny  SMT}. 


Fig.  32.  Spectral  time  history  and 
bathymetry  in  selected  beam  pointing 
directions  (bottom  highlight  here  is  Va- 
vilov  SMT). 
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Fig.  33.  Spectral  time  history  and 
bathymetry  in  selected  beam  pointing 
directions  (bottom  highlight  here  is 
Ostica  Island). 
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Fig.  35.  Single-ping  matched- filtered  returns  over  2r  160  kin;  10  Hz  LFM 
pulse,  (a)  broadside  data;  ( b )  endfire  data. 
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Fig.  36.  Ensemble- aver  aged  ( 1 7  pings)  matched-filtered  returns  over  ^  160  km; 
10  Hz  LFM  pulse.  Lower  curve  is  dB  standard  deviation,  (a)  broadside  data, 
(b)  endfire  data. 
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Fig.  37.  Comparison  of  experimental  cumulative  distribution  function  vs  three 
theoretical  models:  Rayleigh  (curve  R),  log-normal  (curve  LN)  and  Weibull 
(curve  W).  (a)  data  under  test  (stationary  ambient  noise);  (b)  (1-CDF)  vs 
median-normalized  intensities;  log/log  scale. 
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Fig  .1.1  (  -nii/ian.v'ii  ■./  .  sperwieiH.il  .  umuUtiv  distribution  (unction  vs  three 
theoretic*!  mu, lids  Rayleigh  (curve  H  log  norm  *1  ( curve  LS)  and  Weibull 
(curve  IF  I  data  u/l./i-r  lest  (rererhrntlinn  dominated),  (h)  (l-CDF)  vs 

medlAli-n, >riu slued  intensities.  I., g  log  >.  a/o 


Fig.  id  Power  spectra  of  loyarnhmic  intensities,  (a)  ambient  noise  data, 
(b)  reverberation  dominated  data 


Pig.  40.  Comparison  of  experimental  cumulative  distribution  function  vs  three 
theoretical  models:  Rayleigh  (curve  R).  log-normal  (cu.ve  LN)  and  Weibull 
(curve  W).  (a)  data  under  test  (fluctuations  component):  (b)  (l-CDF)  vs 
median -normalized  intensities;  log/log  scale. 
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